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ABSTRACT 

The  geometric  analysis  of  macroscopic  fractures  in  the  Foothills,  Eastern 
Ranges,  and  Main  Ranges  of  the  Canadian  Rocky  Mountains  reveals  a  systematic 
fracure  pattern.  Most  observed  joints  are  found  to  fall  into  either,  tv'O  conjugate 
hkl  sets,  a  Okl  set  replacing  the  hkl  sets,  or  an  hOI  set. 

New  methods  of  fracture  analysis  are  introduced  which  utilize  an  IBM 
1620  computer  system.  The  counting  out  and  plotting  of  orientation  diagrams,  the 
evaluation  of  individual  maxima,  $  ~axis  determinations,  and  the  determination  of 
individual  planar  intersections,  are  all  accomplished  using  computer  programs. 

Principal  stress  trajectories  are  postulated  to  have  been  subparallel  to 
bedding  at  the  time  of  joint  formation.  Changes  in  the  dihedral  angle  between 
conjugate  shear  fractures  reveal  that  areas  presently  at  equal  elevations  were  at 
one  time  subject  to  different  overburden  pressures.  The  joints  are  postulated  to 
have  formed  in  response  to  unloading  by  erosion.  Three  changes  in  the  principal 
stress  directions  during  unloading  account  for  the  orientation  of  the  observed  joint 
sets  and  their  directions  of  movement. 
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INTRODUCTION 


Jointing  has  received  only  a  limited  amount  of  attention  in  structural 
investigations  of  the  Canadian  Rocky  Mountains.  The  analysis  of  mesoscopic 
fractures  has  shed  new  light  onto  the  tectonic  history  of  many  areas,  and  a  study 
of  the  joints  in  the  Canadian  Rocky  Mountains  promised  to  be  a  profitable  project. 
During  the  summer  of  1963,  a  total  of  10,250  joints  were  measured  at  localities 
in  the  Foothills,  Eastern  Ranges,  and  Main  Ranges  of  the  Canadian  Rocky  Mountains, 
The  purpose  of  the  study  may  be  summarized  as  follows: 

a)  to  describe  the  different  types  of  mesoscopic  fractures  found  in  the  area,  to 
separate  these  into  groups  of  apparently  common  origin,  and  to  attempt  to 
determine  the  relative  ages  of  the  groups; 

b)  to  determine  the  regional  preferred  orientation  of  the  most  important  mesoscopic 
fractures; 

c)  to  correlate  the  mesoscopic  geometry  with  the  macroscopic  geometry; 

d)  to  suggest  a  kinematic  and  dynamic  interpretation  of  the  observed  fracture 


pattern . 
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JOINTING  IN  SEDIMENTARY  ROCKS 

Introduction 

By  far  the  most  common  and  yet  least  understood  mesoscopic  planar 
discontinuities  in  both  undisturbed  and  deformed  sedimentary  rocks,  are  the 
fracture  planes  more  or  less  normal  to  bedding,  known  as  joints.  Usage  of  the  term 
joint  varies  somewhat;  before  delving  into  the  various  classifications  of  joints 
and  their  possible  modes  of  origin  the  definition  of  the  term  will  now  be  briefly 
reviewed. 

Billings  (1954,  p.  106),  Hills  (1963,  p.  150),  and  Metz  (1957,  p„  59) 
defined  joints  as  fractures  normal  to  bedding,  or  nearly  so,  along  which  no  visible 
movement  has  taken  place.  However,  since  all  open  joints  have  undergone  move¬ 
ment  at  right  angles  to  their  surface  and  since  joints  are  planes  of  weakness  along 
which  movement  can  readily  occur,  the  definitions  set  forth  by  Turner  (1948,  p.  181), 
Price  (1959,  p.  149)  and  the  A.G,  I.  Glossary  Supplement  (1960,  p,  35)  are  perhaps 
preferable.  The  A.G.  I.  Glossary  Supplement  definition  reads  as  follows:  Joints  - 
fractures  in  rocks,  generally  more  or  less  vertical  or  transverse  to  bedding,  along 
which  no  appreciable  movement  has  taken  place.  " 

Classification 

Geometric  classification.  On  the  basis  of  their  attitude  with  respect  to  bedding, 
joints  have  been  classified  as  strike,  dip,  and  oblique  or  diagonal  joints  (Billings, 
1954,  p.  107).  Similarly  the  terms  longitudinal,  transverse  or  cross  joints,  and 
diagonal  joints  have  been  used  as  geometric  classifications  when  the  basis  of 
classification  was  taken  to  be  the  attitude  with  respect  to  fold  axes  (Hills,  1963, 

pp.  281-282). 
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The  most  satisfactory  geometric  classification  of  joints,  and  the  one  which 
is  followed  in  this  thesis,  is  based  on  a  system  of  orthogonal  fabric  axes  a,  b,  and  c 
which  are  related  to  the  symmetry  elements  on  the  mesoscopic  scale  (Turner  and  Weiss, 
1963,  p.  88).  The  great  flexibility  of  this  classification  and  its  facility  in  the 
correlation  of  minor  and  major  tectonic  elements,  remove  the  ambiguities  of  the 
other  geometric  classifications.  The  designations  of  the  various  types  of  joints 
and  their  equivalents  are  listed  below  (Fig.  1): 

hkO  joints  =  diagonal  joints,  oblique  joints 
hOI  joints  =  strike  joints,  longitudinal  joints 
hOO  joints  =  be  joints 

Ok  I  joints  =  cross  joints,  transverse  joints,  dip  joints 
OkO  joints  =  ac  joints 

hk |  joints  =  joints  not  paralleling  a  fabric  axis 

Genetic  classification.  More  commonly  used  than  the  geometric  classification  is 
the  genetic  classification  of  joints.  As  is  the  pitfall  of  all  genetic  classifications, 
a  mode  of  origin  has  to  be  postulated  before  this  classification  can  be  applied.  A 
great  deal  of  uncertainty  still  exists  as  to  the  exact  origin  of  joints  in  sedimentary 
rocks  and  the  widespread  use  of  a  genetic  terminology  at  the  geometric  stage  of 
structural  analysis  is  unfortunate.  The  two  genetic  classes  of  joints  generally  re¬ 
cognized  are  shear  (hkO)  and  tension  (h0|  or  Ok  I )  joints. 

Shear  joints  are  usually  thought  to  form  obliquely  to  the  principal  stress 
axis.  Commonly  they  are  smooth,  extensive,  and  planar.  They  generally  occur 
in  two  complementary  sets,  one  of  which  may  show  right-hand  displacement  and 
the  other  left-hand  displacement. 

Tension  joints  are  thought  to  form  normal  to  the  direction  of  least  principal 
stress.  They  are  commonly  rough,  warped,  and  show  great  variability  in  strike. 


"  *  ,  •  b  "  ’  1  -■  •«  "•  -  ^  -Jtnioj  ?o  noitcifepl,  „  ,h 
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FIGURE  I 
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The  least  principal  stress  in  the  formation  of  this  type  of  fracture  does  not  necessarily 
have  to  be  tensile  (Griggs  and  Handin,  I960,  pp.  350-351),  but  may  be  com- 
pressional.  The  use  of  the  term  extension  joint  may,  therefore,  be  preferable. 

Morphology 

Joints  which  have  not  undergone  relative  movement  and  have  not  been 
weathered  extensively  often  display  a  number  of  distinctive  surface  features.  The 
classification  and  genetic  significance  of  the  various  structures  found  on  joint 
surfaces  has  recently  been  reviewed  by  Hodgson  (1961c)  and  Robertson  (1961). 

Structures  falling  on  the  main  joint  surface  can  be  differentiated  from  these 
on  the  periphery  or  fringe  of  the  joint.  The  main  joint  face  is  transversed  by  a  central 
axis  from  which  rays  or  plumes  diverge  toward  the  edge  of  the  joint  plane  (Robertson 
1961,  p.  483;  Hodgson,  1961c  p.  494).  The  plumes  consist  of  a  series  of  shallow 
rises  and  hollows  on  the  joint  surface.  The  central  axis  approximately  parallels  the 
upper  and  lower  boundaries  of  the  rock  unit  in  which  it  occurs  (Hodgson,  1961c, 
p.  499).  In  the  peripheral  regions  of  the  joint  plane,  large  numbers  of  small 
en  echelon  fractures  classified  as  f-joints  (Hodgson,  1961c,  p.  501)  or  b-planes 
(Robertson,  1961,  p.  483)  can  be  observed.  Perpendicular  to  the  b-planes  are  a 
series  of  rough,  irregular  fractures  known  as  c-planes  (Robertson,  1961,  p.  484; 
Hodgson,  1961c,  p.  496). 

In  the  opinion  of  some  authors,  the  presence  of  plumose  structures  is 
limited  to  joints  of  "shear"  origin.  (Parker,  1942,  p.  397;  Robertson,  1961,  p. 

486).  Muehlberger  (1961,  p.  215),  on  the  other  hand,  cited  plumose  markings  as 
criteria  for  the  recognition  of  "extension"  fractures. 
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Origin 

Syntectonic  theory  of  origin.  The  theory  of  joint  formation  most  widely  accepted 
among  structural  geologists  postulates  that  joints  in  folded  and  faulted  rocks  can 
be  directly  attributed  to  syntectonic  stresses.  The  existence  of  a  gradual  transition 
from  joints  showing  no  motion  to  faults  is  also  presumed  (de  Sitter,  1956,  p.  122). 
The  orientation  of  the  possible  joint  sets  should,  therefore,  be  predictable  by  the 
Coulomb-Navier  theory  of  fracture.  Complementary  sets  of  shear  joints  should 
develop  parallel  to  the  intermediate  stress  axis  at  an  angle  of  45  -  <f>/ 2  degrees  to 
the  maximum  principal  stress  axis;  tension  joints  perpendicular  to  the  minimum 
principal  stress  axis.  This  indeed  appears  to  be  the  relationship  quite  frequently 
observed. 

Considerable  difficulty  arises,  however,  when  the  time  of  development 
of  the  fractures  is  to  be  determined.  Hoeppener  (1953)  in  his  studies  in  the 
Rheinischen  Schiefergebirge  concluded  that  the  joints  formed  at  the  early  stage 
of  folding.  The  joints  are  found  to  be  displaced  by  bedding-plane  slip.  Hoeppener 
(ibid.)  also  thought  that  these  pre-folding  or  early  -folding  joints  had  been  rotated 
with  respect  to  the  bedding  by  slippage  along  innumerable  microscopic  fractures 
parallel  to  bedding.  Maximum  rotation  was  observed  for  hOO  joints  and  none  for 
OkO  joints.  The  amount  of  rotation  for  given  bedding  dips  can  be  calculated 
theoretically.  According  to  Koebel  (1940),  the  theoretical  values  of  rotation  are 
in  excess  of  those  observed;  Hoeppener  (ibid. )  found  that  the  calculated  and 
measured  values  agreed  reasonably  well  for  bedding  dips  up  to  20  degrees. 

In  the  Swiss  Jura,  Nabholz  (1956)  found  that  hOO  joints  in  one  tight 
anticline  appear  to  have  been  rotated  9  degrees  with  respect  to  bedding,  whereas 
in  three  similar  folds  no  such  rotation  could  be  detected.  The  amount  of  rotation, 
in  this  case,  appeared  to  be  independent  of  the  dip  of  the  beds.  Hoeppener' s 
explanation  for  the  apparent  rotation  of  joints  therefore  does  not  seem  to  hold  in 


'J&91I  :^i;io  r\  ffo:i  -j-'irl 
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all  cases. 

The  variation  in  the  acute  angle  between  sets  of  shear  fractures  is  of  con¬ 
siderable  interest  in  the  interpretation  of  fracture  patterns.  The  Coulomb-Navier 
theory  of  fracture  predicts  that  the  dihedral  angle  between  conjugate  failure  planes, 
regardless  of  the  values  of  the  principal  stresses,  should  be  45  -cj/ 2  degrees, 
for  most  rocks  can  be  approximated  to  be  30  degrees .  One  vould,  therefore,  expect 
the  acute  angle  between  shear  sets  to  be  60  degrees.  Significant  deviations  from 
this  value  have  been  found  by  some  investigators. 

Parker  (1942,  p.  393)  studied  the  jointing  in  central  and  northern  New  York 
and  northern  Pennsylvania  and  found  the  acute  angle  between  two  shear  sets  to  vary 
from  13-30  degrees.  Wager  (1931,  pp.  392-420)  investigated  the  joint  pattern 
in  the  vicinity  of  the  North  Craven  fault  and  observed  that  two  shear  sets  were 
nearly  perpendicular  to  each  other.  In  a  very  gently  folded  anticline  at  Robin 
Hood's  Bay,  Zwart  1951,  p.  3)  found  sets  of  shear  joints  at  an  angle  of  15  degrees 
to  each  other.  Hoeppener  (1953),  in  his  investigation  of  joints,  noted  that  the 
shear  joints  made  an  angle  of  less  than  45  degrees  to  the  perpendicular  of  the  fold 
axis.  Duschatko  (1953,  p.  35)  investigated  the  areal  fracture  pattern  of  the  Lucero 
Uplift  in  New  Mexico  and  found  that  a  single  joint  set  became  replaced  by  two 
conjugate  joint  sets  at  low  dihedral  angle.  The  progressive  increase  in  the  acute 
angle  up  to  65  degrees  was  recorded.  Joint  sets  with  the  largest  dihedral  angles 
were  found  to  be  in  the  structurally  lowest  positions. 

According  to  de  Sitter  (1956,  p.  130),  the  deviation  of  the  acute  angle 
from  the  expected  value  may  be  due  to  "differences  in  rock  properties  or  to  the 
weight  of  the  overburden;  but  on  the  other  hand  it  may  merely  indicate  that  two 
quite  different  kinds  of  shear  joints  exist."  Wager  (1931,  p.  405)  suggested  that  the 
explanation  for  the  variability  may  lie  in  experiments  done  on  natural  materials. 

T.  von  Karman  (191 1)  had  shown  that  the  angle  of  shear  may  be  altered  by  super¬ 
posing  a  hydrostatic  pressure  on  the  non-uniform  pressure  required  to  produce  shear 
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Muehlberger  (1 961 ,  pp,  211-219)  discussed  the  mechanics  of  conjugate 
joint  sets  of  small  dihedral  angel  in  view  of  the  Mohr  theory  of  fracture.  He 
suggested  that  for  low  stresses  the  Mohr  envelope  assumes  a  parabolic  configuration 
and  intersects  the  normal  stress  axis  at  right  angles.  Experimental  data  tends  to 
support  this  assumption.  Under  those  circumstances  there  will  be  a  stress  circle 
with  radius  of  curvature  equal  to  that  of  the  parabola  at  the  vertex.  This  is  the 
largest  stress  circle  for  tensile  failure.  A  slight  increase  in  the  differential  stress 
in  the  aforementioned  case  will  move  the  stress  circle  off  the  vertex  and  it  will 
then  be  tangent  to  the  envelope  at  two  points.  The  resulting  conjugate  fractures 
would  enclose  a  small  dihedral  angle.  Within  a  small  range  of  increasing  stress 
differences,  the  angle  would  increase  to  the  usual  angle  of  shear  failure. 

Post-tectonic  theory  of  origin.  One  of  the  theories  that  attempts  to  explain  the 
formation  of  joints  in  both  undeformed  and  deformed  rocks  postulates  that  unloading, 
or  the  decrease  in  geostatic  pressure,  is  one  of  the  major  factors  in  the  formation 
of  joints.  It  was  noted  by  Hills  (1963,  p.  157)  that,  unless  compensated  for  by 
the  increase  of  density  due  to  compaction,  the  formation  of  tension  joints  involves 
an  increase  in  volume  of  the  affected  rock  mass.  This  increase  ir  volume  suggests 
that  many  joints  may  form  with  decreasing,  rather  than  increasing,  stress  due  to 
relief  of  loading  by  erosion.  Price  (1959,  p.  158)  demonstrated  that  tensile 
stresses,  equal  to  half  the  change  in  gravitational  load,  will  develop  during 
unloading. 

Price  (1959,  pp.  149-167)  used  the  Mohr  theory  of  fracture  to  develop 
a  model  of  events  during  unloading.  The  model  is  based  on  the  assumption  that 
competent  rocks  act  as  brittle  materials  in  the  upper  crust  and  obey  the  elastic 
theory.  Price  (1959,  pp.  156-157)  postulated  that  the  elastic  strain  energy  remains 
stored  in  the  rock  and  residual  stresses  may  exist  which  in  quantity  and  direction 
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are  equal  to  the  stresses  which  acted  at  the  end  of  the  main  tectonic  phase.  Rocks 
are  elastico-viscous  bodies  and  the  residual  stresses  will  tend  to  relax  with  time. 

Very  little  is  known,  however,  about  the  viscosity  of  rocks.  It  may,  therefore, 
be  quite  possible  that  the  relaxation  time  may  be  in  the  order  of  10  million  years. 
This  means  that  the  residual  stresses  would  relax  to  1  %  of  their  original  value 
in  50  million  years  (Price,  1959,  p.  157).  Rocks  may  also  have  a  fundamental 
strength,  that  is,  they  will  not  deform  continuously  unless  a  critical  differential 
stress  is  reached.  Undissipated  stresses  would  then  be  stored  as  residual  stresses. 

Price  (ibid. )  further  showed  that  if  the  initial  conditions  of  the  rock  mass 
approximate  the  hydrostatic  state  the  rocks  actually  go  into  tension  and  develop 
vertical  tension  joints,  provided  uplift  accompanies  erosion.  These  tension  joints 
will  be  randomly  orientated  if  all  horizontal  stresses  are  equal  at  each  depth  level. 
When  horizontal  stresses  are  present,  the  tension  joints  will  form  at  right  angles 
to  the  axis  of  least  principal  stress.  As  soon  as  fracture  has  occurred,  the  least 
principal  stress  is  replaced  by  a  compressive  stress,  and  the  intermediate  principal 
stress  becomes  the  least  principal  stress.  A  second  set  of  tension  fractures  will 
develop  upon  further  uplift.  The  two  sets  will  form  an  orthogonal  system.  !n  cases 
where  the  tectonic  stresses  were  not  relaxed  after  an  initial  period  in  which  the 
greatest  principal  stress  was  horizontal,  the  vertical  load  would  change  from  being 
the  least  principal  stress  to  being  the  intermediate  principal  stress  at  some  level. 

With  further  uplift,  a  set  of  conjugate  shear  fractures  will  form  which  reduces 
the  residual  stresses  considerably.  The  conditions  discussed  in  the  previous 
example  will  be  satisfied  as  uplift  continues  and  the  rocks  may  pass  into  tension 
and  develop  a  system  of  orthogonal  tension  fractures.  In  cases  where  the  residual 
stresses  have  not  been  affected  by  subsequent  tectonic  events,  the  joints  will  still 
faithfully  reflect  the  stress  conditions  as  they  existed  during  the  main  tectonic  phase. 
Equally  plausible  is  the  situation  where  the  stress  conditions  may  become  upset  by 
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subsequent  events,  in  which  case  the  orientation  of  the  joints  bears  no  relationship 
to  the  geometry  of  the  structures  developed  during  the  main  tectonic  phase. 

Formation  of  joints  by  earth  tides.  The  presence  of  well  developed  joint  sets  in 
relatively  undeformed  and  even  flat-lying  sedimentary  rocks  has  led  some  authors 
to  believe  that  the  tension  and  shear  hypotheses  are  inadequate  in  explaining  the 
formation  of  joints.  Hodgson  (1961,  pp.  1-38),  investigating  the  jointing  in  the 
Comb  Ridge-Navajo  Mountain  region,  found  as  many  as  six  well  developed  joint 
sets  in  very  gently  folded  rocks.  The  fractures  bear  no  obvious  geometrical  relationship 
to  the  folding  and  a  number  of  different  stress  systems  would  have  to  be  invoked,  if 
one  were  to  explain  these  fractures  by  compressive  stresses.  Hodgson  (ibid.)  concluded 
that  the  fractures  represent  fatigue  phenomena  resulting  from  the  application  of  external 
stresses.  Most  other  investigators  of  fracture  patterns  in  flat-lying  rocks  also  advocate 
the  formation  of  joints  in  response  of  an  external  stress  field  (Blanchet,  1957,  p.  1754; 
Mol  lard,  1958,  p.  Ill;  Haman,  1961,  p.  9). 

Kendall  and  Briggs  (1933)  were  the  first  to  suggest  that  semi-diurnal  earth 
tides  produce  joints  by  a  torsion  mechanism.  In  addition  to  the  semi-diurnal  lunar 
tides,  Blanchet  (1957,  p.  1755)  also  mentioned  the  possible  effect  of  the  lunar  diurnal 
tides,  the  solar  diurnal  tides,  and  the  solar  semi-diurnal  tides,  each  with  their  own 
period  and  wave  length.  The  effects  of  the  body  tides  on  the  continental  crust  are 
as  foil  ows: 

a)  They  cause  periodic  increases  and  decreases  in  the  porosity  of  near-surface  rocks, 
which  results  in  a  corresponding  rise  and  fall  of  the  local  water  table  (Robinson, 
1939,  pp.  656-666). 

b)  Ihey  cause  the  cyclic  variation  in  the  flow  of  some  hot  springs  (Vorster,  1956). 

c)  They  cause  the  detectable  independent  tilting  of  static  fault  blocks  (Nishimura, 


1950,  p.  359). 
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d)  They  cause  crustal  rise  and  fall  of  a  magnitude  which  has  been  variously  estimated 

at  6  -  21  inches  (Holmes,  1963,  p.  1  41 1 )  and  9  -  1  4  inches  (Blanchet,  1957, 

p.  1754). 

Under  this  hypothesis  joints  are  considered  to  represent  fatigue  phenomena 
resulting  from  the  cyclic  repetition  of  small  stresses.  Earth  tides  act  approximately 
four  times  a  day  or  1  1/2  billion  times  per  million  years.  The  similarity  between 
the  orientation  of  plumose  structures  found  on  joint  surfaces  and  structures  on  the 
surfaces  of  fatigue  fractures  in  metals  and  concrete  was  taken  as  evidence  by  Hodgson 
(1961 ,  p.  23)  for  this  mode  of  origin  of  joints. 

Kendall  and  Briggs  (1933)  postulated  that  earth  tides  controlled  both  the 
formation  of  the  joints  and  their  direction.  The  chief  directions  of  jointing  in  the 
northern  hemisphere  should,  therefore,  be  northeast  and  northwest,  in  coincidence 
with  directions  of  maximum  shearing  stress  caused  by  the  tidal  forces.  Joint  patterns 
do  not  always  show  such  a  pattern  and  modern  workers  have  adopted  the  hypothesis  that 
the  direction,  possibly  even  the  original  formation  of  the  joints,  could  be  regionally 
controlled.  Holmes  (1963,  p.  1411)  pointed  out  that  the  orientation  of  the  developing 
systems  of  fractures  is  probably  controlled  by  independent  stresses  and  that  the  con¬ 
tinuous  tidal  stresses  would  extend  the  fracture  pattern  under  a  minimum  of  internal 
stress.  Hodgson  (1961,  p.  36)  considered  the  possibility  that  the  fracture  pattern 
could  be  inherited  by  each  newly  deposited  rock  unit  from  the  jointed  rock  beneath, 
advocating  thereby  the  upward  migration  of  joints.  In  the  Bright  Angel  area,  Arizona, 
Hodgson  (1961b,  pp.  95-97)  found  that  joints  in  sedimentary  rocks  overlying 
Precambrian  metamorphic  rocks,  in  part,  reflect  the  major  structural  directions  found 
in  the  Archean  metamorphic  complex.  This  seems  to  indicate  that  the  Precambrian 
fracture  pattern  has  been  imposed  on  each  successive  younger  rock  unit. 


12 


COLLECTION  OF  FIELD  DATA 


Selection  of  Localities 


Previous  joint  studies  have  shown  that  the  structural  setting,  lithology, 
and  bedding  thickness  play  important  roles  in  both  the  development  and  frequency 
of  jointing  (Harris,  Taylor,  and  Walper,  1960,  pp.  1859-1860;  Fitzgerald,  1963, 
p,  131).  One  of  the  primary  goals  of  this  investigation  was  to  establish  the  relationship 
between  the  regional  structure  and  the  macroscopic  fractures.  It  was  desirable  to  keep 
both  lithology  and  bedding  thickness  constant.  Other  factors  influencing  the  choice  of 
sampling  locations  included  accessibility,  nature  of  the  outcrop,  and  the  presence  of 
well  developed  fractures. 

The  sandstone  units  in  the  Cardium  formation  of  the  Foothills  of  southwestern 
Alberta  most  nearly  satisfied  these  requisites.  Sandstones  in  the  Cardium  formation 
have  a  fairly  constant  lithology  for  considerable  distances  along  regional  strike, 
although  they  may  change  rapidly  across  it.  Joint  orientations  in  these  sandstone 
units  were  obtained  in  six  anticlines  along  the  Bow  River,  one  anticline  on  the  Red  Deer 
River,  the  limb  of  a  syncline  at  South  Creek,  and  a  monocline  at  Cripple  Creek 

(Fig.  3). 


Work  done  in  the  Eastern  and  Main  Ranges  of  the  Rocky  Mountains  was  on 
a  reconnaissance  basis  to  investigate  the  feasibility  of  more  extensive  studies  in  the 
future.  Joints  were  measured  in  the  competent  units  of  the  Palliser,  Alexo,  Banff, 
and  Rundle  formations  in  a  number  of  different  thrust  sheets  of  the  Eastern  Ranges 
(Fig.  3).  Sampling  in  the  Main  Ranges  was  confined  to  one  locality  situated  in  the 
Eldon  formation  on  the  east  limb  of  the  Bow  River  anticline. 
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Sampling  Procedures 

In  addition  to  the  strike  and  dip,  the  size  and  relative  regularity  of  each 
joint  were  noted.  Also  recorded  were  the  attitude  of  the  slickensides  both  along 
the  joints  and  the  bedding  planes;  the  relative  movement  along  slickensides  wherever 
discernible;  the  attitude  of  the  bedding  planes;  and  the  thickness,  lithology,  and 
stratigraphic  position  of  the  fractured  unit. 

In  order  to  avoid  subjective  selection  in  the  fractures  measured,  all  the 
fractures  which  were  planar,  or  could  be  approximated  by  a  plane,  were  measured. 
Exceptions  were  fractures  smaller  than  a  few  square  inches  in  surface  area,  fractures 
in  the  immediate  proximity  of  known  faults,  curviplanar  fractures,  and  fractures 
definitely  due  to  blasting.  Errors  in  the  measurement  of  the  strike  and  dip  are 
estimated  to  be  between  one  and  five  degrees,  depending  on  the  irregulatiry  of  the 
fracture  surface. 

No  quantitative  measurements  of  the  fracture  frequency  were  made  in 
the  course  of  the  study.  In  order  to  yield  meaningful  results,  joint  frequency  studies 
should  be  carried  out  in  areas  with  large  outcrops  of  bedding  planes  or  at  least  two 
dimensions  of  the  rock  body.  In  the  writer's  opinion  the  exposure  in  the  areas  investi¬ 
gated  did  not  warrant  this  type  of  study.  Only  general  observations  on  the  fracture 
frequency  were  made. 

A  domain,  as  used  in  this  study,  represents  any  three  dimensional  rock 
body  statistically  homogeneous  with  respect  to  the  fracturing.  Outcrops  were  divided 
into  a  number  of  subdomains  of  approximately  10  foot  radius.  Subdomains  were  chosen 
to  have  no  noticeable  changes  in  lithology  or  bedding  attitude.  A  number  of  bedding 
attitudes  and  all  of  the  joints  were  measured  in  each  subdomain.  Subdomains  showing 
no  appreciable  difference  in  the  bedding  or  fracture  attitudes  were  then  grouped  into 
domains  containing  200-300  joint  orientations.  Orientation  diagrams  were  prepared 
for  all  domains. 
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Number  of  Observations  Required 

The  number  of  observations  of  fracture  attitudes  required  to  produce 
valid  orientation  diagrams  poses  a  critical  question  in  any  sampling  program  for 
joints.  Mueller  (1933)  concluded  that  200  observations  were  needed  and  that  possibly 
more  should  be  measured  in  areas  of  unknown  character.  Pincus  (1951,  p.  101) 
considered  80  orientations  per  locality  to  be  the  minimum  value,  while  Spencer 
(1959,  p.  475)  used  1 00-120  measurements.  The  values  quoted  by  Pincus  and 
Spencer  represent  the  minimum  number  of  fracture  attitues  necessary  to  establish  the 
general  pattern,  but  not  the  optimum  number  of  observations  required  for  good  accuracy 
in  the  quantitative  evaluation  of  the  fracture  data. 

In  the  writer's  opinion  the  number  of  observations  required  depends  both 
on  the  number  of  fracture  sets  present  in  a  particular  outcrop,  and  the  dispersion 
shown  by  each  of  these  sets.  The  dependence  of  the  number  of  fractures  needed  in 
the  analysis  on  the  dispersion  can  be  illustrated  by  plotting  the  number  of  fractures 
versus  the  confidence  raius  about  the  mean,  for  a  number  of  fracture  sets  (Fig.  2). 
Fracture  sets  having  different  precision  parameters  (Fisher,  1953,  p.  303),  that  is, 
different  dispersions,  can  be  seen  to  form  independent  plots.  In  the  accompanying 
diagram  set  1  has  a  low  dispersion,  while  sets  2  and  3  have  a  higher  but  equal 
dispersion.  The  plots  show  that  the  number  of  fractures  required  to  reach  a  specified 
value  of  the  confidence  radius  about  the  mean  depends  on  the  dispersion.  It  should 
also  be  noted  that  any  increase  in  the  number  of  observations  will  result  in  an  increase 
in  the  confidence  in  the  mean  fracture  orientations.  An  optimum  number  of  observa¬ 
tions  has  only  a  minimal  effect  on  the  confidence  radius.  In  the  example  illustrated 
(Fig.  2),  this  number  appears  to  be  55  observations  for  set  1  and  65  observations  for 
both  sets  2  and  3.  As  approximately  10  per  cent  of  the  observed  fractures  in  this  area 
do  not  fall  within  one  of  the  fracture  sets,  205  joint  attitudes  should  be  obtained  at 
this  local ity . 


NUMBER  OF  FRACTURE  ATTITUDES 
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Because  the  dispersion  was  one  of  the  parameters  to  be  investigated, 
it  was  found  necessary  to  determine  how  many  observations  per  joint  set  would  be 
required  to  arrive  at  a  consistent  value  of  the  precision  parameter  (Fisher,  1953, 
p.  303).  A  plot  of  the  precision  parameter  of  various  fracture  sets  versus  the  number 
of  fracture  orientations  reveals  that  70  observations  per  set  are  required  to  arrive 
at  a  consistent  value  of  the  precision  parameter.  When  the  dispersion  of  the  fractures 
is  to  be  investigated,  230  fracture  orientations  should  be  measured  at  this  locality. 

The  dispersion  of  the  fractures  can  vary  considerably  from  outcrop  to 
outcrop,  not  all  sets  are  equally  developed  in  any  one  outcrop,  and  the  number  of 
fracture  sets  present  is  not  always  known  beforehand.  It  was,  therefore,  attempted 
to  measure  250-300  joints  per  domain. 


■ 


18 


DATA  PROCESSING 

Introduction 

In  view  of  the  large  number  of  bits  of  information  to  be  processed,  it 
was  found  that  conventional  methods  of  analysis  were  unsatisfactory.  A  number 
of  computer  programs  were  designed  to  aid  in  the  preparation  and  evaluation  of  the 
orientation  diagrams. 

The  equipment  used  consisted  of  the  IBM  1620  Data  Processing  System 
at  the  University  of  Alberta,  Edmonton,  with  1622  Card  Read  Punch.  1623  Core 
Store  Unit  with  one  60K  module  and  a  407  Line  Printer. 

The  programs  were  written  in  FORTRAN  II  and  compiled  with  a  floating¬ 
point  word  length  of  8  characters  and  a  fixed-point  word  length  of  4  characters. 

Program  and  monitor  statements  conform  with  specifications  set  forth  in  the  "FORTRAN 
II  for  the  Magnetic  Tape  1620"  manual  of  the  Computing  Center,  University  of  Alberta, 
Edmonton . 

Preparation  of  Pole  Density  Diagrams:  Program  913107-001 

Introduction .  Joint  studies  in  structural  geology  originated  in  the  latter  half  of  the 
nineteenth  century  in  regions  where  joints  showed  a  great  degree  of  regularity  and 
could  easily  be  defined  by  a  single  observation  of  their  attitude.  Soon  investigators 
found,  however,  that  in  areas,  where  the  joint  pattern  is  less  clearly  defined,  the 
joints  also  exhibit  definite  preferred  orientations.  In  such  cases,  geologists  found  it 
necessary  to  estimate  the  mean  orientation  of  the  joints  in  the  outcrop;  a  method 
which  is  still  widely  used  today. 

Salomon  (1911)  and  his  students  conducted  extensive  joint  surveys  along 
the  Rhine  graben.  By  careful  and  detailed  cartography  they  were  able  to  define 
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regional  patterns  in  the  jointing.  Although  they  took  only  a  few  readings  per  outcrop, 
this  school  laid  the  foundation  for  the  statistical  treatment  of  joints, 

Stiny  (1925,  p.  873)  postulated  that  the  mean  attitudes  of  joints  could  not 
be  estimated  accurately  in  the  field,  and  he  was  the  first  to  advocate  the  measurement 
of  a  large  number  of  joints  per  outcrop,  in  order  to  determine  the  true  joint  direction. 
His  ideas  were  readily  accepted.  Soon  various  methods  were  designed  to  evaluate 
statistical  numbers  of  joint  attitudes.  Strike  frequency  diagrams  and  rectangular 
coordinates  were  frequently  used.  These  methods  were  found  to  be  unsuitable  in 
most  cases-  however,  because  they  only  illustrate  the  variation  in  one  variable, 
whereas  the  attitude  of  any  plane  is  defined  by  two  dependent  variables. 

The  attitude  of  a  plane  in  space  can  be  accurately  represented  on  a  stereo¬ 
graphic  projection.  A  stereographic  net  to  plot  joints  was  first  used  by  Seitz 
(1914/17).  As  equal  areas  on  the  reference  sphere  do  not  remain  equal  on  the 
stereographic  projection,  the  pole  densities  are  distorted  by  this  method.  To  over¬ 
come  this  difficulty,  Schmidt  (1925)  introduced  the  use  of  the  equal-area  projection 
to  petrofabrics.  This  method  of  investigation  was  first  used  in  studying  joints  by 
Rueger  (1928).  In  order  to  facilitate  an  easier  interpretation  of  the  diagrams,  it 
became  customary,  to  count  out  the  diagrams  with  circular  counters  and  to  draw 
density  contours  (Knopf  and  Ingerson,  1938,  pp.  245-252). 

Strand  (1944,  p.  1  12)  pointed  out  that  accurate  contouring  required  a 
change  of  the  unit  counting  area  from  circular  at  the  center  of  the  diagram  to 
elliptical  towards  the  margin.  Duschatko  (1955,  p.  1522)  designed  a  counter  con¬ 
sisting  of  a  number  of  ellipses  with  increasing  ellipticity  toward  the  periphery  of 
the  projection,  in  order  to  eliminate  some  of  the  distortion  introduced  during  the 
counting  out  process.  A  large  number  of  ellipses  would  be  necessary  to  obtain  true 
pole  densities. 


' 

.  4  • 


20 


Statement  of  problem.  Ideally  petrofabric  diagrams  should  not  be  counted  out  in 
their  projection,  but  on  the  reference  sphere.  Pole  densities  obtained  on  the 
reference  sphere  can  then  be  projected  onto  the  equatorial  plane  and  will  faithfully 
reflect  the  true  pole  densities.  Until  recent  years,  this  type  of  analysis  would  have 
been  impossible  for  large  numbers  of  observations.  Problems  involving  repetitive 
intricate  calculations  are  ideally  suited  for  modern  high-speed  digital  computers. 

A  computer  program  was  designed  to  count  out  petrofabric  diagrams  on  the  reference 
sphere  and  to  print  out  the  projections  of  the  densities  onto  the  equatorial  plane. 

Method.  In  order  that  a  number  (NCDS)  of  attitudes  of  planes  or  lineations  can  be 
processed  by  this  program,  each  attitude  has  to  be  in  one  of  the  following  forms: 
strike  (NSTR),  dip  (NDIP),  and  direction  (NDIR);  trend  (NTR)  and  plunge  (NPL); 
or  direction  cosines  (A,  B,  G).  When  in  one  of  the  above  forms,  the  attitude  will 
be  referred  to  as  a  data  reading.  The  computer  converts  the  strike,  dip,  and 
direction  of  a  plane  to  the  trend  (TR)  and  plunge  (PL)  of  the  normal  to  the  plane. 

All  the  trends  and  plunges  are  then  changed  into  direction  cosines  A  yjLA  .,  and  y. 


of  unit  length  using  the  relationships: 

A  .  =  cos  TR.  *  cos 

i  i 

PL 

i 

a) 

fA.  =  sin  TR.  *  cos 

PL. 

i 

(2) 

Y\=sin  PL. 

(3) 

The  counting  out  of  the  poles  in  the  reference  sphere  is  accomplished  by 
the  use  of  a  circular  cone  with  semivertical  angle  of  o(.  The  value  otoC  is  determined 
by  the  value  of  P,  the  percentage  of  the  total  area  of  the  projection  which  the 
counter  is  required  to  cover.  The  value  of  P  and  the  number  (N)  of  counting 
locations  have  to  be  chosen  such  that  the  counting  cones  overlap  in  such  a  manner 
as  to  assure  that  all  points  on  the  sphere  are  counted.  At  the  same  time  P  should 
be  kept  small  in  order  to  obtain  the  best  possible  definition  of  any  differences  in 
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the  densities  of  the  points.  In  most  petrofabric  studies  P  has  been  taken  to  be  1 .0 
and  N  as  31 7. 

The  centre  of  the  counting  cone  is  placed  at  specified  and  previously 
calculated  locations.  If  the  direction  cosines  of  the  counting  locations  are 
\  c'  A  c'  anc*  J/ c'  anc^  readings  X  y  y  ar|d  then  the  angle 

between  their  directional  lines  Ss  given  by 

ft  <4> 

If  ((X  the  data  reading  fall  s  within  the  counting  cone.  The  number  of  data 

readings  in  each  data  set  falling  within  the  counting  cone  at  each  counting  locality 
is  noted,  and  converted  into  the  percentage  of  the  total  number  of  points  processed. 

In  order  to  facilitate  a  simple  output  format  using  the  IBM  407  Line  Printer, 
the  counting  locations  were  chosen  to  form  a  regular  grid  with  a  spacing  of  R/10 
(R  =  radius  of  the  reference  sphere)  on  the  equatorial  plane.  For  the  convenience 
of  accurate  contouring,  these  were  augmented  by  1  4  additional  points  near  the 
periphery  of  the  projection.  The  counting  points  are  projected  only  from  the  lower 
hemisphere  of  the  reference  sphere  onto  the  equatorial  plane  by  equal-area  projection. 


Limitations  on  parameters. 

NCDS  ^  500 
0  ^  NSTR  ^  180 

0  ^  NDIP  90 

0  ^  NDIR  <  2 

0  ^  NTR  ^  360 

0  ^  NPL  ^  90 

-1.0  ^  A  =5:  1.0 

-1.0  ^  6  *--=•  1.0 

0.0  ^  G  1.0 


' 
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Input,  (a)  Parameter  card 


Column 

Content 

Format 

1  ™ ■  m  m  . 

“ - -  “  ■ 

1 

blank 

IX 

i 

CM 

NCDS 

13 

5 

blank 

IX 

6 

Code 

11 

7 

blank 

IX 

o 

i 

00 

Run  No. 

13 

The  following  code  has  been  employed  in 

column  6  of  the  parameter  card 

0  if  the  input  data 

consists  of  the  strike. 

dip,  and  direction  planes; 

1  if  the  input  data  consists  of  the  trend  and  plunge  of  1  inears; 

2  if  the  input  data  consists  of  the  direction  cosines  of  lines  (output  from 

Program  913107-003  or  004). 

(b)  Data  cards 

Data  cards  differ  in  format  according  to  the  code  used  in  column  6  of  the 

parameter  card  as  follows: 

Code  =  0 

Column 

Content 

Format 

1 

blank 

IX 

2-4 

NSTR 

13 

5 

blank 

IX 

6-7 

NDIP 

12 

8 

blank 

IX 

9 

NDIR 

11 

10  -  15 

blank 

7X 

16  -  80 

can  be  used  for 

coding  and  sorting. 

- 
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The  direction  of  dip  (NDIRj  of  planes  has  been  coded  in  the  follov'ing 

manner: 

0  if  the  dip  is  90  degrees; 

1  if  the  dip  direction  is  N,  NE,  NW,  or  W; 

2  if  the  dip  direction  is  S,  SE,  SW,  or  E. 

N  and  S  are  only  valid  for  strikes  of  90  degrees,  and  similarly  E  and  W  for  strikes 
of  0  and  180  degrees;  otherwise  the  NE,  NW,  SE,  or  SW  quadrants  have  to  be 
specified.  For  example,  179/85  W  is  an  invalid  data  reading  and  should  read 
1 79/85  SW. 


Code  =  1 

Column 

Content 

Format 

1 

blank 

IX 

2-4 

NTR 

13 

5 

blank 

IX 

6-7 

N  PL 

12 

8-15 

blank 

10X 

16  -  80 

can  be  used  for 

coding  and  sorting. 

Code  =  2 

Column 

Content 

Format 

1  -  7 

blank 

7X 

8-18 

A 

F11.8 

19-29 

B 

F11.8 

Co 

o 

1 

6 

G 

F11.8 
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(c)  Termination  card 

Content  Format 

blank  1 4X 

5  11 


Column 
1  -  14 
15 


Order  of  cards. 

(a)  Program  913107-001  (pp.  25-27) 

(b)  Direction  cosines  of  counting  locations  (Appendix  I) 

(c)  Parameter  card 

(d)  Data  cards 

Repeat  (c)  and  (d)  for  each  subsequent  run. 

(e)  Termination  card 

Output.  The  output  of  this  program  consists  of  two  grids  of  numbers,  one  representing 
the  per  cent  density  of  poles  per  counting  cone,  and  the  other  the  actual  number  of 
points  counted  per  counting  cone.  The  grids,  when  circumscribed  by  a  circle  of 
10  inch  diameter,  represents  the  densities  of  points  at  counting  points  on  the 
reference  sphere  that  have  been  projected  onto  the  equatorial  plane. 

The  numerals  are  right-hand  orientated  and  contours  should  be  placed 
through  the  center  of  the  right-hand  digit. 

Test  data  and  results.  The  accompanying  output  sample  (pp.  28-29)  represents 
237  attitudes  of  joint  planes  measured  on  the  east  limb  of  Kananaskis  anticline. 

The  contoured  output  (Fig,  5)  of  the  computer  can  be  compared  with  a  similar 
diagram  (Fig,  4)  prepared  by  the  conventional  graphic  methods  using  a  circular 
counter.  In  general  the  agreement  is  excellent  in  the  central  regions  of  the  diagram. 
Marked  differences  exist  in  the  peripheral  regions  if  the  graphical  counter  used  is 
circular  rather  than  elliptical. 
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..I  9??!07  POINT  nFNSlTY  DIAGRAMS  G.K.MUECKE 

..I  PLEASE  SAVE  CARDS. 

..load  FORTRAN  execute  DUMP 

DIMENSION  COSA( 333 ) *COSR( 333 ) *COSG( 333 ) *ND( 333) 
DIMENSION  A ( 500 ) * B ( 500 ) *G ( 500 ) 

1  READ  INPUT  TAPF  5*50»NDENS 

DO  2  1  =  1*  NDENS 

2  READ  INPUT  TAPE  5  *  5 1  *  COS A ( I ) > COSB (  I )  ♦ COSG ( I  ) 

READ  INPUT  TAPE  5*52>KPCA 

D  A  =  K  P  C  A 

TEST= 1 . -DA* .01 

THET  =  ATANF( SORTF ( 1./ ( TEST*T EST )-l.  )  ) 

3  DO  4  1  =  1*  NDENS 

4  ND (  I  )  =  0 
NZ  =  0 

P I =3 • 1 41 5 926  . 

R 180  =  1. /I  80. 

CF  =  P I *R 1 8  0 
I  X  =  1 

READ  INPUT  TAPE  5  *  5 0  * NCDS * NZ *  NS TN > NUT S 
I E ( NUTS )45  *46*45 

45  CALL  EXIT 

46  FNCDS=NCDS 
RNCDS=1 00. /FNCDS 
NZ  =NZ+1 

GO  TO  ( 47  *47  *  1200 )  *NZ 

47  DO  12  1=1  *  NCDS 
GO  TO  ( 5  *  10  )  ,NZ 

5  READ  INPUT  TAPE  5  *  5 3 > NSTR , ND I P * ND I R 
PL=90-NDI P 

I F ( NDIR  —  1  )  8  *  7  *  8 

7  TR=NSTR+90 
GO  TO  11 

8  T  R  =NS  T  R-9  0 

I  F  (  T  P  )  9  ,  l  i  ,  1 1 

9  TR=TR+360. 

GO  TO  11 

10  p E AD  INPUT  TAPE  5>53*NTD*NPL 
Tp=NTR 

PL=NPL 

11  PL=PL*CF 
TR  =  T  R*CF 

A ( I )=COSF(TR)*COSF(PL) 

B (  I  )  =S  I  NF ( TR )*COSF (PL ) 

G ( I ) =SINF ( PL ) 

12  CONTINUE 
GO  TO  1202 

1  200  DO  1201  I  =  1  *  NCOS 

1201  READ  INPUT  TAPE  5  *  5 1  *  A (  I  ) * B (  I ) * G (  I ) 

1202  P M A X  =  0 . 

DO  124  I  =  1 ♦ N  C  D  S 
I F ( PM A X-G (  I )  )  123  > 124* 1 24 

123  PM AX=G ( I ) 

124  CONTINUF 
PMAX=PMAX+TFST 
DO  15  1=1  *  NDENS 


. 

. 
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CSG=COSG( I ) 

I F ( PMAX-CSG ) 15*15,125 
125  CSA=C0SA(I) 

CSB=COSB( I ) 

DENS=0 • 

AT=CSG+TEST 

bt=csg-test 

DO  14  J=1 ,  NCOS 
GE  E  =  G ( J ) 

I F ( GEE-BT ) 14*122*121 

121  IF (AT-GEE)14*122*122 

122  PEST=ABSF( CSA*A( J ) +CSB*8 ( J ) +CSG*GEE ) 
IF (TFST-PEST ) 13*13*14 

13  PFNS=DENS+1. 

14  CONTINUE 

15  NO (  I  ) =DENS 


'<  L  I  K  = 

1 

vVB  I  T  E 

OUT  PUT 

tape 

6 

*  5  8  *  N  D  (  1  ) 

WR  I  TE 

OUT  PUT 

tape 

6 

*66*ND( 318 ) 

*  ND ( 3 1 9 )  *  ND ( 320 ) ,ND(321 ) 

WRITE 

OUTPUT 

TAPE 

6 

*68  *NSTN 

WR  I  TE 

OUTPUT 

TAPE 

6 

*  59* ( ND(  I ) , 

1=2*10 ) 

W  RITE 

OUTPUT 

TAPE 

6 

*57 

WR  ITE 

OUTPUT 

TAPE 

6 

*57 

WR  ITE 

OUTPUT 

TAPE 

6 

,60, ( ND( I ) , 

1=11*23) 

WR  ITE 

OUTPUT 

TAPF 

6 

,  57 

WR  I  TF 

OUTPUT 

TAPE 

6 

*  5  7 

WR  I  TF 

OUTPUT 

tape 

6 

*61 , (ND( I ) , 

1=24*38 ) 

WR  I  TE 

OUTPUT 

TAPF 

6 

*57 

W  RITE 

OUTPUT 

tape 

6 

*57 

W  RITE 

OUTPUT 

tape 

6 

*62*  <  ND(  I ) * 

1=39,55 ) 

WR  ITE 

OUTPUT 

TAPE 

6 

,57 

WRITE 

OUTPUT 

TAPE 

6 

,57 

WRITE 

OUTPUT 

TAPE 

6 

,62, <ND( I ) 9 

1=56,72 ) 

WR  ITE 

OUTPUT 

TAPE 

6 

,  57 

WR  ITE 

OUTPUT 

TAPE 

6 

,57 

WR  ITE 

OUTPUT 

TAPF 

6 

,63, ( ND( I ) , 

1=73,91 ) 

WR  I  T  F 

OUTPUT 

TAPF 

6 

*57 

WR  I  TE 

OUTPUT 

TAPF 

6 

*  57 

WR  ITE 

OUTPUT 

tape 

6 

*63, (ND( I ) * 

1  =  92  *110) 

WO  j  TF 

OUTPUT 

TAPE 

6 

*  57 

WR  ITE 

OUTPUT 

TAPE 

6 

*  57 

W  R  I  T  E 

OUTPUT 

tape 

6 

*63, (ND( I ) , 

1=111*129) 

WR  I  TE 

OUTPUT 

TAPE 

6 

,67*ND( 322 ) 

,ND( 323 ) 

WRITE 

OUTPUT 

TAPE 

6 

,57 

WRITE 

OUT  PUT 

tape 

6 

*63* ( NO ( I ) * 

1=130*148 ) 

WR  ITE 

OUTPUT 

TAPE 

6 

*  67 » ND ( 324 ) 

*  ND ( 325 ) 

WR  ITE 

OUTPUT 

TAPE 

6 

*  57 

WR  I  TF 

OUTPUT 

TAPE 

6 

, 64* ( ND( I ) * 

1=149,169) 

WR  ITE 

OUTPUT 

TAPE 

6 

*  57 

WR  I  TF 

OUTPUT 

TAPF 

6 

,67,ND(326) 

,ND( 327 ) 

WR  I  T  F 

OUTPUT 

TAPF 

6 

,63, ( ND( I ) , 

1  =  1-70, 1  88  ) 

WRITE 

OUTPUT 

TAPF 

6 

*  57 

WR  I  T E 

OUTPUT 

TAPE 

6 

*  67  *  ND ( 328 ) 

*  ND ( 329 ) 

WR  I  T  E 

OUTPUT 

tape 

6 

,63* ( ND( I ) * 

1=189,207) 

WR  I  TE 

OUTPUT 

TAPE 

6 

,57 

WR  ITE 

OUTPUT 

TAPF 

6 

,  57 

.  ■  • 

’ 


I  *  (  I  ) ' 


T  \  qw 


17 

18 

1  9 

20 

21 

22 

50 

5  1 

52 

5  ? 

5  4 

57 

58 

59 

60 

6  1 

62 

63 

64 

65 

6  6 

67 

68 

69 

7  0 
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WRITE 

OUTPUT 

TAPE 

6 , 63, ( ND ( I ) , 1=208 , 226 ) 

WRITE 

OUTPUT 

TAPE 

6,57 

WR  I  TE 

OUTPUT 

tape 

6  *  57 

WR  I  TE 

OUTPUT 

TAPE 

6*63, ( ND ( I ) , 1=227,245 ) 

WR  I  T  E 

OUTPUT 

TAPF 

6  *  57 

WRITE 

OUTPUT 

TAPF 

6  *  57 

WR  I  TE 

OUTPUT 

TAPE 

6*62* ( ND  < I ) ,  1  =  246,262) 

WR  I  TE 

OUTPUT 

tape 

6*57 

WR  I  TE 

OUTPUT 

TAPE 

6,57 

WR  I  TE 

OUTPUT 

TAPF 

6*62* ( ND ( I ) , 1=263,279) 

WR  I  TF 

OUTPUT 

TAPF 

6,57 

W  P  I  T  E 

OUTPUT 

TAPF 

6  >  57 

WR  I  T E 

OUTPUT 

TAPF 

6*61 , ( ND ( I )  *  1=280*294) 

WR  I  TE 

OUTPUT 

TAPF 

6 ,57 

WP  I  TF 

OUTPUT 

TAPF 

6  *57 

WR  I  T  E 

OUTPUT 

TAPF 

6*60, ( ND ( I ) , 1=295,307) 

GO  TO 

(17*18) 

*KL  IK 

WR  I  TE 

OUTPUT 

TAPE 

6*54 

GO  TO 

19 

WRITE 

OUTPUT 

TAPE 

6,69 

WRITE 

OUTPUT 

TAPE 

6  *  70 , KPCA 

WR  ITE 

OUTPUT 

TAPF 

6*  50, ( ND(  I )  ,  1=308 , 316 ) 

WR  ITE 

OUTPUT 

TAPE 

6*57 

WR  I  T  E 

OUTPUT 

tape 

6*66*ND(330) *  N  D ( 3  3 1 )  ,ND(332) 

WP  I  T  F 

OUTPUT 

TAPF 

6  *  6  5  *  ND ( 317) 

GO  TO 

( 20,22 ) 

*KL  IK 

DO  21 

1=1  ,333 

XND=ND  (  I  ) 

ND ( I )=XND*RNCDS+.5 
!<  L  I  K  =  2 
GO  TO  16 
GO  TO  3 

FORMAT ( IX, I3*1X,I1*1X, 13*15) 

FORMAT ( 7X  *  3F1 1 • 8  ) 

FORMAT (13) 

FORMAT ( IX , 13* IX, 12  » IX, 1 1 ) 

FORMAT (21H  NUMBER  OF  POINTS) 

FORMAT ( IX ) 

FORMAT ( 21  HI  DENS  I T Y  D I STR I  BUT  I  ON  * 2  9X , I  5 ) 
F  ORMAT ( 3 0 X  ♦  9  I  5 ) 

FORMAT ( 20 X  ,  1 3  I  5  ) 

FORMAT ( 1 5  X  ,  1 5  I  5  ) 

FORMAT ( 1 0 X  ♦  1 7  I  5 ) 

FORMAT ( 5X  ,  19  I  5  ) 

FORMAT (2115) 

FORMAT ( 50 X , I  5 ) 

FORMAT ( 39X, 15, 16* 1 1 0  *  16) 

FORMAT (  16  *93X,  15 ) 

FORMAT ( !2H  STATION  NO. *14) 

FORMAT ( 2 1 H  PERCENTAGE  OF  POINTS) 

FORMAT ( 4H  PER* I3*13H  PERCFNT  AREA) 

END 
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Determination  of  Mean  Fracture  Planess  Program  913107-002 

Introduction ,  The  quantitative  evaluation  of  orientation  diagrams  has  received 
considerable  attention  in  recent  years,,  In  fracture  analysis  the  determination  of 
the  mean  fracture  planes,  which  are  represented  by  the  observed  maxima,  is  of 
utmost  importance.  These  have  previously  been  determined  by  empirically  choosing 
the  peak  value  of  the  concentration,  or  by  estimating  centers  of  gravity  (Pincus, 

1951,  p.  105;  Spencer,  1959,  p.  480),  Where  peak  values  are  not  symmetrically 
centered,  significant  errors  may  be  introduced  by  this  method. 

Statement  of  problem.  Individual  maxima  can  be  analyzed  using  the  statistical 
treatment  of  the  dispersion  of  points  on  a  sphere,  as  has  been  developed  by  R.A. 
Fisher  (1953).  This  method  of  analysis  has  proved  to  be  a  useful  took  in  the  analysis 
of  paleomagnetic  data,  Cox  and  Doe!!  (1960,  p,  669)  suggested  that  this  method 
of  analysis  could  also  be  applied  to  other  problems  involving  sets  of  vectors  or 
lines.  The  mean  value  defined  by  the  maximum,  the  confidence  radius  about  the 
mean,  and  the  dispersion  can  be  calculated. 

Method,  The  data  processed  by  this  program  consists  of  a  number  (NCDS)  of 
fracture  attitudes  specified  by  their  strike  (NSTR),  dip  (NDIP),  and  direction  (NDIR), 
These  attitudes  are  converted  by  the  computer  to  direction  cosines  (  >  y  y  Y  p 
by  using  equations  (1),  (2),  and  (3). 

When  using  Fisher's  statistics  one  assumes  that  the  sample  was  randomly 
drawn  from  a  population  in  which  the  points  are  distributed  with  axial  symmetry 
about  the  true  mean  direction,  and  that  the  density  vector  decreases  with  increasing 
angular  displacement  ^  from  the  mean  according  to  the  probability  density  function 

p  =  (K/47f'sinh  K)  exp  (K  cos^) 


(5) 
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K  is  known  as  the  "precision  parameter"  and  is  a  measure  of  the  dispersion 
of  the  points.  K  is  large  for  a  tight  grouping  about  the  mean,  and  zero  for  a 
uniform  distribution  of  the  points  on  the  sphere. 

Fisher  (1953,  p.  296)  showed  that  the  best  estimate  ( A  , ju ,  /)  of  the 
true  mean  direction  is  the  vector  sum  of  the  individual  directions  (  ^  .). 

That  is 


Saj 

II 

A 

Ki 

li 

ii 

M 

(6) 

R 

R 

R 

r2  =  (JA  ,)2  +  (2fx  ,)2  +  (2y.)2 

(7) 

The  azimuth  and  the  inclination  of  the  resultant  vector  are  determined  using  the 
relationships  (Cox  and  Doell,  1960,  p.  669)s 

tan  ^  i  sin  I  =  $  i  (8) 

R  R 


where  D  is  the  declination  and  I  the  inclination  of  the  resultant  vector, 
r  r 

The  accuracy  of  the  calculated  mean  direction  (  A  ,  y)  can  be  estimated 
by  the  cosine  of  the  angle  between  this  mean  and  the  true  mean.  Fisher  (1953, 
p.  303)  showed  that  at  a  probability  level  of  (1  -  P)  the  true  mean  direction  of  the 
population  lies  within  a  circular  cone  about  the  resultant  vector  with  a  semi- 
vertical  angle  <X ,  given  for  K>  3  by  the  expression 


cos  <X  =  1  -  N  -  R 
R 


(9) 


where  N  is  the  total  number  of  vectors  included  in  the  analysis.  Usually  P  is 
chosen  to  be  .05,  which  means  that  there  is  one  chance  in  twenty  that  the  true 
mean  direction  lies  outside  the  cone  of  confidence. 


. 
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Where  3,  Fisher  (ibid, )  showed  that  "k",  the  best  estimate  of  K,  is 

given  by 

k  -  N  -  1 

N-R  (10) 

As  N  is  increased  without  limit,  this  estimate  of  the  precision  parameter  approaches 
the  true  value  of  K,  o(  on  the  other  hand  becomes  infinitely  small  (Cox  and  Doell, 
I960,  p.  671). 

Wilson  (1959,  p,  751)  pointed  out  that  Fisher's  statistics  are  based  on  an 
assumed  Gaussian  distribution  of  density  of  the  measured  vectors  about  the  mean 
direction.  Flowever,  this  type  of  distribution  is  not  always  valid,  and  in  a  rigorous 
analysis  the  observations  should  first  be  fitted  to  Fisher's  distribution.  In  cases 
where  there  is  no  good  fit,  any  circle  of  confidence  containing  the  mean  direction 
with  a  specified  probability  may  contain  falsely  implied  precision.  Wilson  (1959, 
p„  755)  suggested  the  use  of  the  root  mean  square  deviation,  given  by 

=  (2  (N-R)  /  N  ),/2  (11) 

and  the  standard  deviation  of  the  mean  direction,  given  by 

i  =  ^/N1/2  (12) 

In  some  cases  it  may  be  advantageous  to  use  (f)  and  £  ,  because  no  assumptions 
about  the  distribution  of  the  line  or  vector  densities  have  to  be  made. 

Because  Fisher's  statistics  are  designed  for  vectors,  lines  have  to  be  con¬ 
verted  to  vectors  before  proceeding  with  the  analysis.  This  can  be  accomplished 
by  dividing  the  reference  sphere  into  two  hemispheres  along  the  equatorial  plane, 
and  arbitrarily  assigning  a  positive  sense  to  all  the  lines  falling  into  the  lower 
hemisphere.  In  many  petrofabric  problems  planes  are  steeply  dipping  and  lineations 
gently  plunging.  The  poles  of  the  same  maximum  will  therefore  commonly  fall  at 
diametrically  opposed  positions  on  the  hemisphere.  In  order  to  account  for  this 
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behaviour  in  the  computer  program,  a  dextral  (HIM)  and  a  sinistral  (NUM)  azimuthal 
limit  had  to  be  drawn  for  each  maximum  (Fig.  4). 

Limitations  on  parameters. 


NCDS 
0  <  P 
180  <  LLIM 
0  <-  NUM 
0  ^  NSTR 
0  ^  NDIP 
0  NDIR 

Input,  (a)  Parameter  cards. 

Card  1 
Column 
1 

2-4 

5 

6-7 

8 

9-11 

12 

13 

Card  2 
Column 
1 

2-4 

5 

6-8 


^  500 

'S'  10 

<  360 
^  360 

4  180 

90 


Content  Format 

blank  IX 

Station  No.  13 

blank  IX 

Maximum  No.  12 

blank  IX 

NCDS  13 

blank  IX 

P  11 

Content  Format 

blank  IX 

LLIM  13 

blank  IX 

NLIM 


13 


- 
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(b)  Data  cards 


Column 

Content 

Format 

1 

blank 

IX 

2-4 

NSTR 

13 

5 

blank 

IX 

6-7 

NDIP 

12 

8 

blank 

IX 

9 

NDIR 

11 

Order  of  cards. 

(a)  Program  913107-002  (pp.  36-37) 

(b)  Parameter  card  1 

(c)  Parameter  card  2 

(d)  Data  cards 

Repeat  (b)  to  (d)  for  subsequent  runs.  Each  run  of  data,  except  the  first,  must 
be  preceeded  by  the  following  cards:  . .  I 

..  BRANCH  0040R 

Output.  The  program  puts  out  the  following  information: 

(a)  Station  number  and  number  of  the  maximum  under  consideration. 

(b)  Number  of  data  readings  included  in  the  analysis. 

(c)  Declination  of  the  normal  to  the  mean  fracture  plane  in  degrees  of  azimuth. 

(d)  Inclinat  ion  of  the  normal  to  the  mean  fracture  plane  in  degrees  from  the  horizontal. 

(e)  An  estimate  of  the  precision  parameter. 

(f)  Confidence  radius  in  degrees  about  the  mean  direction,  for  a  specified  probability 
level . 

(g)  Root  mean  square  deviation  in  degrees. 

(h)  Standard  deviation  from  the  mean  in  degrees. 


■  i!  . 


b  r  .jp?  osm  too#  lf ) 

'c  ic  'eb  n  ♦<?  (H) 
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..I  ° 1 3 1 u 7  MEAN  FRACTURE  PLANE  MODIFIED  G  •  K  • MUECKF 
.  .  I  PLEASE  SAVE  DUMP. 

.  .  I  PLEASE  SAVE  CARDS. 

..LOAD  FORTRAN  EXECUTE 

DIMFNS ION  A( 500 ) *R ( 500 ) >G ( 500 ) 

P I =  3 • 1 4T  5°26 
R  1  80= 1  . /180. 

CF=PI*R!80 

READ  50»NSTN»NDEN>NCDS*NP 
P=NP 

READ  51 >LLIM»NLIM 
S A  =  0 . 0 
SB=0 . 0 
SG=0 • 0 

1  DO  12  I =1 *  N  C  D  S 

READ  52*NSTR»NDIP*NDIR 
PL  =  R0-ND I P 
IF (NDJR-1 ) 3,2,3 

2  TR=NSTR+90 
GO  TO  8 

3  TR  =  NSTR  —  90 

I  F  (  T  R  )  4 , 5  >  5 

4  T  R  =  T R  +  360 . 

5  DLIM=LLIM 

IF ( TR-DLIM)6*10*10 

6  SL  I  M  =  NL  I  M 

IF (SLIM-DLIM)8»9,9 

8  IF ( TP-90.  )  10,10,9 

9  TR  =  TR  +  180. 

PL  =  PL*CF 
To  =  TP*CF 

A ( I )=COSF(TR)*COSF(PL) 

•R (  I  ) =S INF ( TP )*COSF ( PL ) 

G( I ) =-S I NF ( PL ) 

GO  TO  11 

10  PL=PL*CF 
TR  =  TR*CF 

A (  I ) =COSF ( TR ) *C  OSF ( PL ) 

R ( I ) =S I NF ( TP )*COSF ( PL ) 

G (  I  ) =S INF ( PL  ) 

11  $  A  =  S  A  +  A (  I  ) 

S  B  =  S  P  +  B  (  I  ) 

12  SG  =  SG+GU) 

R=SQRTf( ( SA*SA)+(SB*SR)+(SG*SG) ) 

SININ=SG/R 
T  ADE  =  SB/S A 

REIN=ATANF(SI N IN/SQRTF ( 1 • -S I N I N*S ININ) ) 

RE IN  =  RE  IN/CF 
RE  IN  =  ABSF (REIN) 

REDE=ATANF( TADE ) 

REDE=REDE/CF 
I F ( SG ) 121 *122*122 

121  S A=- ( S A ) 

S  B  = - ( SR ) 

122  IF(SA)lr',13,]R 

13  IF(SB)16*17,17 


15 

1  6 

1  7 

c;  . 

r,  i 

52 

6  C 

5] 

6  2 

6  3 

6  4 

6  5 

66 

67 


37 


P  E  DE  =  1 8  0 . +RFDE 
GO  TO  17 
REDE=360.+RFDF 
CDS=NCDS 

X  <  = ( COS-1 .  ) / ( COS-R ) 

PUNCH  6  C  >  N  S  T  N  *  N  D  F  N 
P O' N r H  6 NCOS 
PUNCH  A  ^ , P  F  0  F 
PUNCH  6 2  »  RF  I  N 
PUNrH  61  »  XK 
PERC= ( 1 Ou  * -P ) 

P  =  P  * .  0  ] 

Q=l. /(COS-1. ) 

RR= ( CDS-P ) /P 
PP= ( ( 1 . /P ) **0 ) -1 . 

pp  =pp*pp 

COS AL= 1 .-RR 

F  U  N  =  A  Q  S  F (  ( 1 .-rOSAL*COSAL ) / ( COSAL )  ) 

FUN=SOPTF ( FUN ) 

A.  L  =  A  T  A N  P  (  FUN  ) 

A  L  =  A  L  /  r  F 

OFLT A  =  SOR  T  F ( 2 .* ( COS-R ) /COS ) 

F  T  A  =  OF L T A / SOR T F ( CDS ) 

OELTA=OFLTA/CF 
ETA=ETA/CF 
PUNCH  64*PEPC*AL 
PUNCH  66»DELTA 
PUNCH  6  7  »  E  T A 

FORMAT  (  1  X  »  I  3  »  1  X  ♦  I  2  *  1  X  <»  I  3  *  1  X  *  1 1  ) 

FORMAT ( ] X  *  I  3  *  1 X  *  13) 

FORMAT ( a  X  *  I  8,  i X»  I  2  *  1 X,  II  ) 

FORMAT (//PH  S  T  a  T 1 0  N I  4  *  1 4  H  MFAN  F^ACTURFJ? ) 

FORMAT  (  2 3 H  PRrClSION  p/\pamfTFR»12X>F8.3) 
f o R M AT ( 2 2 H  RFSULTANT  I MCL  I  NAT  I  ON ♦ 1  OX » F 8 • 3 ) 

FORMAT ( 22 H  RESULTANT  OFCL I N A  T I  0  N  *  1 0  X  »  F  8 • 3 ) 

F  O  p  M  A  T ( 1 X  > I  2  » 1 X  *  2  P  H  PER  CENT  CONF I DENCF  RAO  I  US  FP.3) 
FORMAT ( // 13H  mFASURFMENTS  123 ) 

FORMAT  (  14H  °‘AS  OEV  I  AT  J  ON  ♦  1  8X  »  F  8 . 3  ) 

FORMAT ( 19H  STANDARD  DFV I  AT  I  ON  *  1 3 X > F3 . 3 ) 

CALL  p X  I T 
END 


- 
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Test  data  and  results.  An  example  solved  by  Fisher  (1935,  p.  304)  which  consists  of 
nine  paleomagnetic  directions  determined  by  Hospers  (1953)  was  chosen  to  test  the 
program.  The  answers  derived  by  Fisher  and  the  computer  are  identical,  (pp.  39-40) 


Determination  of  Beta  Axes  and  the  Preparation  of  Beta  Diagrams: 

Program  913107-003 

Introduction.  The  most  pronounced  tendency  of  any  particular  set  of  s-surfaces  in  a 
tectonite  body  is  to  intersect  along  a  common  axis.  Where  the  s-surfaces  are  folded 
cyl indrical ly ,  these  intersections  tend  to  cluster  around  a  mean  value,  known  as  the 
-axis,  which  defines  the  axis  of  folding  B  (Turner  and  Weiss,  1963,  pp.  154-155). 
Conventional  methods  of  analysis  involve  the  determination  of  lines  of  inter¬ 
sections  on  an  equal-area  net,  and  the  contouring  of  the  diagram.  The  peak  value 
or  center  of  gravity  of  the  resulting  maximum  is  then  taken  to  define  the  -axis. 

Statement  of  problem.  It  is  desired  to  determine  a  large  number  of  planar  intersections, 
to  calculate  their  mean  direction,  and  to  produce  an  equal-area  projection  showing  the 
distribution  of  these  intersections. 


Method.  The  data  processed  by  this  program  consists  of  a  number  (NCDS)  of  attitudes 
of  s-surfaces  specified  by  their  strike  (NSTR),  dip  (NDIP),  and  direction  (NDIR).  These 
attitudes  are  converted  to  direction  cosines  (A  *,  y  j)  using  equations  (1),  (2), 
and  (3). 


The  normal  equation  of  any  plane  passing  through  the  origin,  whose  normal 
has  the  direction  cosines  ^  ,  /f,  and  y ,  is  given  by 

A  x  +^y  +y  z  =0  (13) 


For  two  given  planes 

^  +  X  ]z  =  0 

and  A  2X  +  +  2Z  =  ® 


(14) 

(15) 


the  direction  cosines  of  their  line  of  intersection  are  given  by  A  , /A  ,  and 


■ 
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. .  I 
. .  I 

..I  913107-002 
..BRANCH  0040R 
222  02  009  5 
343  062 
073  24  2 
152  21  2 
127  20  2 
11  7  0  8  2 
089  11  2 
'J  96  17  2 
140  21  2 
088  31  2 
134  39  2 


INPUT  DATA 

TEST  DAM  BY  FISHER. 
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.  .  I  PROGRAM  OUTPUT 

•  •  I 

..I  913107-002  TEST  DATA  BY  FISHER. 
..BRANCH  004CR 


STATION  222  MEAN  FRACTURE  2 


MEASUREMENTS  9 

RESULTANT  DECLINATION  24.269 

RESULTANT  INCLINATION  70.771 

PRECISION  PARAMETER  35.131 

95  PER  CENT  CONFIDENCE  RADIUS  8.755 
RMS  DEVIATION  12.888 

STANDARD  DEVIATION  4.296 
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It  can  then  be  shown  that 


AS  Yi  "A  2  Y  \ 


06) 


=  Z4  l  <3^2  -^2  o] 

'  s  - - - ~ — _ 

K 

(17) 

jA  .=  ^  2  h  -^1 

K 

(18) 

Y  ^  ^  i  2  -  A  2/^1 

0  ^ 

K 

(19) 

where  K  =  {jA  ^  2  /  ,)2  +  (>2^ ,  -  ^  ]  / 2)2  +  {  ^  ^  2  '  ^  2^  / 

(20) 


The  (}  -axis  is  then  found  by  calculating  the  vector  sum  of  all  the  individual 
directions  and  using  equations  (6),  (7),  and  (8).  The  confidence  radius  about  the 
|3”axis  and  the  precision  parameter  are  calculated  from  equations  (9)  and  (10) 
respectively. 

The  number  of  mutual  intersections  of  N  non-parallel  planes  is  given  by 


N(N  -  l)/2 


(21) 


It  should  be  noted  that  the  program  will  put  out  the  direction  cosines  of  any  specified 
number  (IF INI)  of  intersections.  The  format  of  this  output  is  compatible  with  the  input 
specifications  for  Program  913107-001  .  The  output  of  the  present  program  can, 
therefore,  be  used  directly  as  data  cards  in  the  counting  out  of  a  beta  diagram. 


Limitations  on  parameters.  The  limitations  are  the  same  as  for  Program  913107-002, 
with  the  exception  of  0  ^  IF  INI  ^  9999. 


3i  sd  'i'JOJ, 

o  ]A\:  )  >dmun 
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Input,  (a)  Parameter  cards 

Card  1 


Column 

Content 

Format 

1 

blank 

IX 

i 

CM 

Station  No. 

13 

5 

blank 

IX 

6-8 

NCDS 

13 

9 

blank 

IX 

10 

P 

11 

11 

blank 

IX 

12  -  15 

IF  1 N 1 

14 

Card  2  has  the 

same  format  as  in 

Program  913107-002. 

(b)  Data  cards  have  the  same  format  as  in  Program  913107-002. 

Order  of  cards. 

(a)  Program  913107-003  (pp.  43-45) 

(b)  Parameter  card  1 

(c)  Parameter  card  2 

(d)  Data  cards 

Repeat  (b)  to  (d)  for  each  subsequent  run.  Each  run  of  data,  except  the  first, 
must  be  preceded  by  the  following  cards: 

..I 

..BRANCH  0040R 

Output.  The  following  information  is  put  out  by  the  program: 

(a)  direction  cosines  of  a  specified  number  of  intersections 

(b)  station  number, 

(c)  number  of  s-surfaces  used  in  the  analysis. 
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..I  913107  BETA  DETERMINATIONS  G. K • MUECKE 

..I  S A V F  CARDS.  PLFASE  RECORD  RUNNING  TIME 
..LOAD  FORTRAN  EXFCUT E 

DIMENSION  A (300)  »B ( 300 )  >G(300)  *  A  I  NT (  300 ) »BI NT ( 3  00 )  >GINT ( 300 ) 
PI =3. 1415926 

R 1 80  =  ]  .  / 1 8 0 . 

CF=PI*R180 

READ  50*NSTN*NCDS*NP> IFINI 
P=NP 

READ  51  >LLIM,,NLIM 
5  A  =  0 . 0 
SB=0 • 0 
S  G  =  0 . 0 
K=  0 

1  DO  11  I  =  1  *  NCDS 

READ  52>NSTP»NDIP*NDIR 

PL=90-NDIP 

I F ( ND I R-l )  3  *  2  *  3 

2  TR=NSTR+90 
GO  TO  10 

3  T  R  =NSTR-9  0 
IF(TR)4*10*10 

4  TR  =TR  +  360 . 

10  TR=TR*CE 
PL.  =  PL*CE 

A (  I  )  =COSF ( TR )*COSF ( PL ) 

P (  I  )  =S  INF ( TR ) *COSF ( PL ) 

11  G (  I  )  =S  I  NF ( PL ) 

20  DO  21  J  =  1  * NCDS 

AA= A ( J ) 

BA=B ( J ) 

GA=G ( J ) 

GA=G ( J ) 

3  0  DO  3 R  L  =  1  »  NCDS 
A  B  =  A  (  L  ) 

9  B  =  R  (  L  ) 

GB=G( L ) 

I F ( AA-AR ) 33  »  31 >  33 

31  IF ( RA-BP) 33*32,33 

32  IF (GA-GP) 33*37*33 

33  DTRE=(AA*BB)-(AB*BA ) 

DONE= ( B  A  *  G  B ) -( EB*GA ) 

IF (DTRE)37*332*331 

332  I F ( DONE ) 37*331 >331 
331  DTWO=( AR*GA)-(AA*GB) 

D I  V  I  =  SQRT E (  ( DONE*DONF )  +  ( DTWO*DTWO)  + ( DTRE*DTRE)  ) 

K  =  K  + 1 

A  I  NT ( L ) =DONE/D I VI 
R I  NT ( L ) =  D T WO / D I V  I 
GI NT ( L ) =DTRF/DI VI 
IE  ( K  —  I R I N I ) 40*40*41 

40  PUNCH  71 * AINT ( L )  *B  INT ( L )  *GINT ( L ) 

41  I F ( A  I  N  T ( L )  )33p*337>33° 

338  TRINT=ATANF(BINT(L)/AINT(L) ) /C F 
I F ( A  I  NT ( L )  ) 335*333,333 

333  IF ( BINT (L )) 334*336*336 
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334  TR  INT=360.+TRINT 
GO  TO  336 

^  3  ^  TP INT=1 80.+TRINT 
GO  TO  336 

33^  I F ( BINT ( L )) 339*340*340 

339  TR INT=90. 

GO  TO  336 

340  TR INT=270. 

336  DLIM=LLIM 

IF ( TRINT-DLIM) 34*38*38 
34  S  L  I M  =  N  L  I M 

IF (SLIM-DLIM)36*36*36 

36  IF ( T  R I N  T- 9  0 .  )3 8 *38*36 
3  6  A  I  NT ( L ) =-4INT ( L ) 

R I  NT (L ) =-B  I  NT ( L  ) 

GI NT ( L ) =-GINT ( L ) 

GO  TO  38 

37  AINT(L)=0.0 
BINT ( L  )  =0.0 
GI NT ( L  )  =0.0 

38  SA  =  SA  +  A  INT ( L ) 

SB=SB+B INT ( L ) 

39  SG  =  5G+GINT(L  ) 

21  A  A  =  A ( J ) 

R  =  SORT  F ( ( SA*S  A )  +  ( SB*SB )  +  ( SG*SG )  ) 

S I NIN=SG/R 
T  ADE  =  SB/S A 

RF ! N  =  aT ANF ( S I N IN/ SORT F ( 1 • - S I N I N # S I N I N  )  ) 
RE IN=REIN/CF 
RE IN=ABSF ( REIN ) 

R  E 0 E  =  A  T A, N F  (  TADE  ) 

REDE=REOE/CF 
IF(SG) 121*122*122 

121  S  A  =  - ( S  A ) 

S B  =  - ( SB  ) 

122  IF ( SA) 1 6* 1 4* 14 
14  IF (SB) 17* 18* 18 

1 6  pEDF=l RO.+RFDE 
GO  TO  18 

17  REDF  =  36  0 . +REDE 

18  Z I N  T  =  K 
XK=(ZINT-1. ) / (ZINT-P ) 

PUNCH  6 0 , NST N 

PUNCH  66* NCOS 
PUNCH  66  *  K 
PUNCH  6 1  *  X K 
PUNCH  6 3  *  REDE 
PUNCH  62* RE  IN 
P=  P*. 01 

0= 1 . / ( ZINT-I .  ) 

RR= ( Z INT-P ) /R 
pp= (  ( 1 . /P ) **Q) -1  . 

RO  =RR*PP 

C  0  vS  A  L  =  1  •  -  R  R 

F  U  N  =  A  B S F (  (  1  . -COS AL*COSAL ) / ( COSAL )  ) 
FUN=SORTF ( FUN ) 


':t£  ! 


50 

51 

52 

60 

61 

62 

63 

6  4 

65 

66 

71 
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AL  = ATANF ( FUN ) 

AL= AL/CF 
PUNCH  64 *  AL 

FORMAT ( IX*  I3>1X» I3»1X. II *1X*I4) 

FORMAT (  IX*  I3*1X» 13 ) 

FORMAT { IX* I3*1X* I2*1X* II ) 

FORMAT ( // 8 H  STAT I0NI4*4X* 24H  BETA  AXIS  DETERMINATION) 
FORMAT ( 2 1 H  PRECISION  PARAMETER  F10.3) 

FORMAT ( 23H  RESULTANT  INCLINATION  F8.3) 

FORMAT ( 23H  RESULTANT  DECLINATION  F8.3) 

FORMAT ( 19H  CONFIDENCE  RADIUS  F12.3) 

FORMAT ( // 1 3H  MEASUREMENTS  1 14) 

FORMAT ( 1 4H  I  NT EPSECT I ONS 1 1 3 ) 

FORMAT ( 7X  *  3F 1 1 • 8 ) 

CALL  EXIT 
END 


■  . 

•  .  •  ' 
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(d)  number  of  intersections  calculated  and  used  in  calculating  the  mean, 

(e)  precision  parameter, 

(f)  azimuth  of  the  delination  of  the  j]  -axis  in  degrees, 

(g)  inclination  of  the  (3  -axis  in  degrees  from  the  horizontal, 

(h)  confidence  radius  about  the  (3  -axis  for  a  specified  probablility  level. 

Test  data  and  results.  Four  planes  which  intersect  along  a  common  line  on  the 
stereographic  net  were  chosen  as  test  data.  The  attitude  of  their  intersection  agrees 
closely  with  the  value  calculated  by  the  computer,  (pp.  47-48) 

Determination  of  Cleavage- Bedding  Intersections:  Program  913107-004 

Statement  of  problem.  The  intersection  between  two  different  sets  of  s-surfaces  may 
sometimes  yield  useful  structural  information.  The  intersection  between  bedding  and 
axial  plane  cleavage  may  be  used  to  define  the  axis  of  folding  B,  or  it  may  be  desired 
to  find  the  lines  of  intersection  of  different  fracture  sets. 

Method.  The  data  processed  by  this  program  consists  of  a  number  (NCDS)  of  attitudes 
of  s-surfaces  specified  by  their  strike  (NSTR),  dip  (NDIP),  and  direction  (NDIR). 
These  attitudes  are  converted  to  direction  cosines  (X  yJA  y  Y us'n9  equations 
(1),  (2),  and  (3). 

Intersections  are  calculated  using  equations  (17),  (18),  (19),  and  (20).  The 
program  differs  from  913107-003  in  that  intersections  of  chosen  planes  are  calculated, 
rather  than  all  the  possible  intersections  of  a  set  of  s-surfaces. 

Limitations  on  parameters. 

NCDS  300 
0  NSTR  180 

0  NDIP  90 


0 


NDIR 


2 


I 
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. .  I 
. .  I 

..I  913107  BETA 
•.BRANCH  0040R 
999  0 C 4  5 
1  R  0  2  7  0 
170  30  2 
146  26  2 
086  46  2 
054  90  0 


INPUT  DATA 
DETERMINAT IONS 


G.K.MUECKE 
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..I  PROGRAM  OUTPUT 

•  •  I 

•  •I  9 1 ^ 1 0 7  BFTA  DETrRMI NAT  IONS  G.K*MUECKE 

..BRANCH  0040R 


STATION  999  BETA  AXIS  DETERMINATION 


MEASUREMENTS  4 

INTERSECTIONS  6 

PRECISION  PARAMETER  1510.026 
RESULTANT  DECLINATION  235.695 
RESULTANT  INCLINATION  26.388 
CONFIDENCE  RADIUS  1.724 
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Input,  (a)  Parameter  card 

Column  Content  Format 

1  blank  IX 

2-4  Station  N  |3 

5  blank  IX 

6-8  NCDS  13 


(b)  Data  cards  have  the  same  format  as  in  Program  913107-002.  Each  data 
set  should  consist  of  a  series  of  pairs  of  data  cards,  each  pair  being  a 
bedding  attitude  followed  by  a  cleavage  attitude. 

Order  of  cards. 

(a)  Program  913107-004  (pp.  50-51) 

(b)  Parameter  card 

(c)  Data  cards 

Repeat  (b)  and  (c)  for  subsequent  runs.  Each  run  of  data,  except  the  first, 
must  be  preceded  by  the  following  cards: 

..I 

..BRANCH  0040R 

Output.  The  following  information  is  put  out  by  the  program: 

(a)  station  number, 

(b)  intersection  number, 

(c)  bedding  attitude, 

(d)  cleavage  attitude, 

(e)  trend  of  the  cleavage-bedding  intersection, 

(f)  plunge  of  the  cleavage-bedding  intersection, 

(g)  direction  cosines  of  the  cleavage-bedding  intersection. 

It  should  be  noted  that  the  program  puts  out  the  direction  cosines  of  the  inter- 


o  liocj  i  >r>9  dob  ' .  i 
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..I  913107  BEDDING-CLEAVAGE  INTERSECTIONS  G.K.MUFCKE 
•.I  PLEASE  SAVE  CARDS. 

..LOAD  FORTRAN  DUMP 

DIMENSION  A ( 300 ) »  p ( 300 ) >G(300) »  KSTp ( 700 ) »  KD  I  P ( 300 ) >KD I R ( 300  ) 
PI =3. 1415926 
R 1 80=1 . /I 80 . 

C  F  =  P I *R 18  0 
READ  1 00  *  NSTN  »  NCDS 
PUNCH  200  , NSTN 
I F IN I=NCDS-1 
K  =  0 

1  DO  11  1=1, NCDS 

READ  101 ,NSTR,NDIP,NDIR 
KSTR ( I ) =NSTo 
KD I P ( I ) =ND I P 
KD IP ( I ) =ND I R 
PL=90-NDI P 
IF (NDIR-1 ) 3 , 2  *  3 

2  TR=NSTR+90 
GO  TO  10 

3  TR=NSTR-90 

IF ( TR)4,10,10 

4  TR=TR+360. 

10  TR=TR*CF 
PL=PL*CF 

A (  I ) =COSF ( TP ) *COSF ( PL  ) 

B (  I  ) =S INF ( TP )*GOSF ( PL  ) 

11  G( I )=SINF( PL) 

20  DO  21  J  =  1 *  I F I N I »  2 
A  A  =  A  (  J  ) 

BA=B( J ) 

G A  =  G ( J ) 

JSTR=KSTR ( J ) 

JD  I  P  =  KD I P ( J ) 

JD  IP  =  KDIR ( J ) 

L  =  J  +  l 
AB=A ( L ) 

BB=B( L ) 

G  B  =  G  (  L  ) 

LSTR=KSTR ( L ) 

LD  I  P  =  KD I P ( L ) 

L  D  I  R  =  KD I R ( L ) 

I F ( AA-AB ) 33, 31 ,33 

31  IF(BA-BB)33, 32*33 

32  I F ( GA-GB ) 33,43 ,33 

33  DTRE= ( AA*BB )-( AB*BA  ) 

DONE= ( BA*GR ) - ( B  B  *  G  A ) 

D  T  W  0  = ( AB*GA ) - ( AA*GB ) 

D I  V  I =SQRT  F (  ( DONF*DONE )  +  ( DTWO*DTWO )  +  ( DTRE*DTRE )  ) 

K  =  K+l 

A  I mT  =  DONF /D I  V I 
PI NT=DTWO/DI VI 
GI NT=DTRF/DI VI 
T  I  P  =  A  I N  T 
T  AP=BI NT 
T OE  =  G I  NT 


.Odf+^  T  ST  C» 
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IF ( AINT  ) 36*43*36 
36  TRINT=ATANF(BlNT/AINT)/CF 
GO  TO  30 

40  I F ( B  I  N T  )  4 1  *  4 2  *  4 2 

41  T  p i n  T  =  q  0 . 

GO  TO  420 

42  TPINT=270, 

42  0  PL  INT  =  ATANF( GINT /SQRTM 1 . -G I NT*G I NT )  ) 

PL  INT  =  ABSF( PLINT/CF ) 

GO  TO  65 

43  A  I N  T  =  0 • 0 
R I NT  =  0  •  0 
GI NT=0.0 

5  0  PL  IMT  =  ATANF(GIMT/SQPTF( 1 .-GINT*GTNT )  ) 

PLINTsABSFtPLINT/CF) 

I  F  (  G I  M  T  )  6  0  *  6  1  *  6 1 

6  0  ai,NT  =  ~AINT 

R I NT  =  -B I  NT 

61  I F ( A  I  NT )  6  3  >  62  *  62 

62  IF ( BINT )64*65*65 

63  TRINT=180,+TRINT 
GO  TO  65 

64  TR  INT  =  360.  +  TRINT 

65  PUNCH  2 0 1  *  K 

PUNCH  202 > JSTR* JDIP* JDIR 
PUNCH  203*LST5>,LDIP>LniR 
PUNCH  2G4»TPINT 
PUNCH  205  »  P L I N T 
21  PUNCH  206  * T I P , TAP  *  TOE 
100  FORMAT ( IX* I  3* IX* 13 ) 

1  C  1  FORMAT  (  ].X*I3>1X*  I2*1X*  II  ) 

2  Ou  FORMAT ( // 8  H  ST AT  I  ON  I  4  * 4X , 30H  CLEAVAGF-BEDO I NG  INTERSECTION) 

231  FORMAT ( // 1 7 H  INTERSECTION  NO. *13) 

232  FORMAT ( / 1 7  H  BEDDING  ATT  I TUDE * 7X » I  3  *  1 X > I  2  *  1 X  *  I  1 ) 

203  FORMAT (18H  CLEAVAGE  ATT  I TUDE * 6X , I  3 » 1 X  *  I  2  *  IX , 1 1 ) 

204  FORMAT ( 20H  INTERSECTION  TREND  F11.3) 

205  FORMAT ( 2 1 H  INTERSECTION  PLUNGE  F10.3) 

206  FORMAT ( 7X  ,3F1 1 .8 ) 

CALL  EXIT 

END 


r  >'■  *  ^  ! 
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sections  in  a  format  compatible  to  the  input  for  Program  913107-001,  The  direction 
cosines  can,  therefore,  be  used  directly  to  produce  orientation  diagrams. 

Test  data  and  results,  A  data  set  consisting  of  four  pairs  of  bedding  and  cleavage 
attitudes  respectively  was  used  to  test  the  program.  Results  obtained  by  the  computer 
program  check  with  those  obtained  using  the  stereographic  net.  (pp.  53-54) 

Preparation  of  Uneation  Diagrams 

A  method,  which  permits  the  illustration  of  I  inears  and  their  relationship 
to  the  plane  in  which  they  occur,  has  been  devised  by  Hoeppener  (1957,  pp.  19-20). 
The  attitude  of  the  si ickensides,  the  joint  plane  in  which  they  occur,  and  the 
relative  movement  can  be  shown  in  one  diagram. 

The  attitude  of  any  plane  passing  through  the  center  of  the  reference 
sphere  can  be  defined  by  the  penetration  pole  of  the  normal  to  the  plane.  A  linear 
lying  in  the  plane  can  be  defined  by  a  family  of  great  circles  in  which  the  linear  lies. 
There  is,  however,  only  one  great  circle  which  is  common  to  both  the  penetration 
pole  of  the  plane  and  the  penetration  pole  of  the  linear.  This  unique  great  circle 
plus  the  penetration  pole  of  the  plane  clearly  show  the  relationship  between  the  plane 
and  the  linear.  On  the  projection  it  is  usually  sufficient  to  give  the  great  circle 
only  in  the  vicinity  of  the  penetration  pole  on  the  sphere.  The  sense  of  relative 
movement,  if  known,  can  then  be  given  by  marking  the  direction  of  movement  of  the 
hanging  wall  by  an  arrow. 
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•  I 

INPUT 

•  I 

.1  913107  TEST 

data 

.BRANCH  004QR 

999  008 

.160  16  2 

RDG 

1  0  8  2  0  1 

CLE 

020  78  2 

BDG 

064  80  1 

CLE 

055  24  2 

BDG 

135  28  2 

CLE 

170  36  2 

BDG 

070  50  2 

CLE 

DATA 
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..I  PROGRAM  OUTPUT 

.  .  I 

..I  P 1 3 1 07  bfdding-cleavagf  intfpsfgt IONS  G.K 


STATION  9QQ  CLEAVAGE-BFODING  INTFPSECTION 


I  INTERSECT  ion  NO.  1 


BEDDING  ATTITUDE  160  16  2 


CLEAVAGE  ATTITUDF  108  20  1 

INTERSECTION  TREND  310.687 

INTERSECTION  PLUNGE  7.991 

.64559932  -.75091602  .13902055 


INTERSECTION  NO.  2 

BEDDING  ATTITUDE  20  78  2 

CLEAVAGE  ATTITUDE  64  80  1 

INTERSECTION  TREND  44.155 

INTERSECTION  PLUNGE  62.551 

-.33070879  -.32110202  -.88742617 


INTFPSECTION  NO.  3 

PENDING  ATTITUDE 
CLEAVAGE  ATTITUDE 
INTERSECTION  trend 
INTERSECTION  PLUNGE 
-•93795415 


INTERSECT  ION  NO.  4 
P FDD  I NG  ATTITUDE 

cleavage  attitude 
intersection  trend 
INTERSECTION  plunge 
.65104118 


55  24  2 
135  28  2 
178.977 
20.264 

>1673619  .34635526 


170  36  2 
70  30  2 
221.502 
29.623 

7603689  -.49429434 


.MUECKE 
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GENERAL  GEOLOGY 

Introduction 

Joints  were  measured  in  three  of  the  structural  sub-provinces  of  the  Canadian 
Rocky  Mountains  (Fig.  3).  The  Foothills  sub-province,  bound  to  the  east  by  the 
undisturbed  Interior  Plains  and  to  the  west  by  the  McConnell  fault,  is  composed 
mainly  of  Mesozoic  sedimentary  rocks.  Closely  spaced  west-dipping  thrust  faults 
and  tight  folds  are  the  most  prominent  macroscopic  structural  elements.  The  Eastern 
Ranges,  which  border  the  Foothills  to  the  west,  consist  predominantly  of  Upper 
Paleozoic  strata.  Imbricate  thrust  sheets,  showing  usually  only  a  minor  amount  of 
folding,  predominate  in  this  region.  The  Castle  Mountain  thrust  forms  the  eastern 
boundary  of  the  Main  Range  sub-province.  Rocks  in  the  Castle  Mountain  thrust 
sheet  are  of  Lower  Paleozoic  age  and  form  gently  undulating  structures. 

Bow  River  Area 

Previous  work.  The  general  geology  along  the  Bow  River  between  Kananaskis  Falls 
and  Morley  was  previously  studied  by  Cairnes  (1914),  Rutherford  (1927),  and 
Beach  (1943).  A  limited  number  of  observations  on  the  mesocscopic  structures  were 
made  by  Roeder  (1960,  p.  591)  in  the  vicinity  of  Kananaskis  Falls. 

The  writer's  investigations  were  confined  to  a  seven  mile  long  strip  along 
the  Bow  River,  east  of  Kananaskis  Falls  (Fig.  6). 

Stratigraphy.  All  beds  observed  in  the  course  of  this  study  belong  to  the  Upper 
Cretaceous  Alberta  group,  which  can  be  divided  into  the  following  formations: 

Upper  Cretaceous  Wapiab?  formation  1775' 

Cardium  formation  220’  -  230 


Blackstone  formation 


700’  -  800 


. 


FIGURE  6 
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Fracture  measurements  were  carried  out  in  the  Cardium  formation,  which  can 
be  further  subdivided  using  distinctive  changes  in  lithology.  A  section  of  the  Cardium 
formation  measured  by  Cairnes  (1914,  pp.  27-28)  at  the  mouth  of  Oldfort  Creek  may 
be  summarized  as  follows: 

Overlying  beds  -  Blackstone  formation 
Cardium  formation 


Upper  sandstone 

unit 

70 

Dark  shale  unit 

60 

Middle  sandstone  unit 

30 

Dark  shale  unit 

20 

Lower  sandstone 

unit 

40 

Underlying  beds  -  Wapiabi  formation 

The  middle  sandstone  unit  thins  to  the  east  and  is  represented  only  by  an 
argillaceous  sandstone  in  the  eastern  portions  of  the  region. 

Lithology  of  sandstone  units.  The  sandstones  of  the  Cardium  formation  are  brownish 
grey  to  grey,  rusty  red  to  brown  weathering,  generally  fine  grained,  quartzose,  and 
well  cemented.  Cross-stratification  in  all  the  sandstone  units  was  found  to  be  more 
common  than  previously  reported.  Bedding  thicknesses  are  1/2  inch  to  1  foot  in  the 
upper  sandstone  unit,  1  foot  to  4  feet  in  the  middle  sandstone  unit,  and  1/2  inch  to 
3  feet  in  the  lower  sandstone  unit. 

Structural  setting.  The  area  «s  located  in  the  severely  disturbed  Foothills  sub-province. 
The  resistant  sandstone  members  of  the  Cardium  formation  produce  readily  traceable 
ridges  north  and  south  of  the  Bow  River,  and  excellent  outcrops  are  found  along  the 
banks  of  the  Bow  River.  The  beds  are  deformed  into  folds  which  trend  north  30  degrees 
west  and  plunge  to  the  northwest.  Six  of  the  anticlines  were  investigated  in  detail. 
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Three  thrust  faults,  striking  from  north  30  degrees  west  to  north-south,  can  be  inferred 
from  the  abrupt  terminations  of  Cardium  outcrops.  Thrust  planes,  although  not 
exposed,  are  thought  to  dip  steeply  to  the  southwest. 

Red  Deer  River  Area 

The  Red  Deer  anticline  is  situated  north  of  the  junction  of  the  Red  Deer 
River  and  the  Alberta  Forestry  Trunk  Road  (Fig.  7).  The  area  lies  within  the 
Foothills  and  its  stratigraphy  is  identical  to  that  discussed  for  the  Bow  River  region. 
Joints  were  measured  in  the  upper  and  lower  sandstone  units  of  the  Cardium  formation. 
25  feet  of  the  upper  and  20  feet  of  the  lower  sandstone  unit  are  exposed.  The  sand¬ 
stones  are  grey  in  color,  rusty  weathering,  fine  grained,  well  sorted,  and  quartzitic. 
Cross-stratification  is  common.  The  thickness  of  the  beds  varies  from  2  inches  to 
2  1/2  feet  in  the  upper  unit,  and  from  6  inches  to  3  feet  in  the  lower  unit. 

The  beds  are  in  a  large  anticlinal  fold  which  has  a  well  exposed  eastern 
limb  and  only  poor  exposures  in  the  western  limb.  The  fold  is  asymmetric  and  its 
axial  plane  dips  steeply  to  the  southwest.  The  anticlinal  axis  trends  approximately 
north  40  degrees  west  and  plunges  10  degrees  to  the  northwest. 

Cripple  Creek  Area 

Erdman  (1946)  has  investigated  the  general  geology  of  this  section  of  the 
Foothills.  Fractures  were  measured  in  the  Cardium  Formation,  which  is  375  feet 
thick  at  this  locality.  The  Cardium  formation  can  be  subdivided  into  the  following 
units  (Erdman,  1946,  p.  2): 

Overlying  beds  -  Wapiabi  formation 
Cardium  formation 


Upper  sandstone  unit 

75 

Middle  sandy  shale  unit 

200 

Lower  sandstone  unit 

100 

Underlying  beds  -  Blackstone  formation 
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The  sandstones  are  grey,  rusty  weathering,  fine  grained,  and  quartzitic.  Cross- 
stratification  and  ripple  marks  are  present,  and  the  bedding  thickness  ranges  from 
4  inches  to  2  1/2  feet. 

The  sampling  locality  is  situated  in  a  southwest-dipping  monoclinal  thrust 
sheet.  The  thrust  shows  a  maximum  stratigraphic  throw  os  5000'  a  few  miles  south 
of  the  area  investigated  (Erdman,  1946,  p.  3). 

South  Creek  Area 

Joints  were  measured  in  the  lower  sandstone  unit  of  the  Cardium  formation. 
The  sandstones  attain  a  thickness  of  60'  at  this  location  and  are  brownish  grey  in 
color,  rusty  brown  to  orange  brown  weathering,  fine  grained,  and  show  frequent 
cross-stratification.  Beds  range  from  2  inches  to  4  feet  in  thickness. 

Fractures  were  measured  at  two  locations  two  miles  apart,  but  in  identical 
structural  positions  on  the  east  limb  of  a  syncline  trending  north  35  degrees  west. 

Exshaw  Area 

Fracturing  in  the  Exshaw  thrust  sheet  of  the  Eastern  Ranges  was  investigated 
at  one  locality  on  Highway  No.  1A  (Fig.  8).  The  resistant,  massive  limestones  of 
the  upper  part  of  the  Palliser  formation  are  exposed  along  the  road  and  show  excellent 
fracture  development.  The  lithology  is  a  dark  grey,  light  grey  weathering  bio- 
pelmicrites. 

The  Exshaw  fault  in  this  area  dips  approximately  60  degrees  to  the  west 
(Clark,  1954,  p.  44).  Beds  in  the  thrust  sheet  show  a  monoclinal  dip  of  40  to  65 
degrees  to  the  west. 


V 


FIGURE  8 
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The  Gap  Area 

The  upper  portion  of  the  Banff  formation  in  the  Lac  des  Arcs  thrust  sheet 
(Fig.  8)  shows  some  fracturing  in  a  road  outcrop  along  Highway  No.  1A.  The  beds 
are  mainly  dark  grey,  medium  grained,  argillaceous  biocalcarenites,  and  dark  grey 
calcilutites.  Frequent  stringers  of  black  chert,  up  to  3  inches  thick,  are  found  near 
the  top  of  the  exposure.  Heavy  calcite  veining  sub-parallel  to  the  bedding  is 
common. 

The  Lac  des  Arcs  thrust  sheet  of  the  Eastern  Ranges  has  been  folded  along 
axes  trending  nearly  parallel  to  the  trace  of  the  fault.  Fracture  orientations  were 
obtained  on  the  western  limb  of  an  anticline  (Fig.  8). 

Roche  Miette  Area 

The  area  investigated  belongs  to  the  Disaster  Point  thrust  sheet  of  the 
Eastern  Ranges.  This  area  has  previously  been  mapped  by  Mountjoy  (1959).  Fractures 
were  measured  in  the  massive  beds  of  the  Devonian  Palliser  formation.  The  lithology 
consists  of  dark  grey,  crypto-  to  fine  crystalline,  dolomitic  limestones.  Bedding 
thickness  varies  from  6  inches  to  5  feet. 

The  beds  lie  in  the  overturned  west  limb  of  an  anticline.  Although  the 
Disaster  Point  fault  dips  to  the  east  in  the  vicinity  of  the  sampling  location,  it  is 
asssumed  to  dip  to  the  southwest  at  depth. 

Morro  Peak  Area 

Road  outcrops  along  the  Jasper-Edmonton  highway  in  the  Colin  thrust  sheet 
show  extensive  fracturing.  Separate  sets  of  measurements  were  made  for  the  Palliser 
and  the  Alexo  formation.  The  Palliser  formation  is  composed  of  dark  grey,  crypto- 
to  fine  crystalline,  slightly  argillaceous  limestones.  Bedding  ranges  from  1  inch  to 
4  feet  in  thickness.  The  Alexo  formation  consists  predominantly  of  yellowish  grey. 


63 


calcareous  siltstones  In  beds  from  6  inches  to  5  feet  thick.  The  beds  of  the  Colin 
thrust  sheet  show  a  uniform  dip  to  the  southwest  in  this  area. 

Mount  Murchison  Area 

The  Mount  Murchison  area  is  situated  in  the  Main  Ranges  of  the  Rocky 
Mountains  and  forms  part  of  the  Castle  Mountain  thrust  sheet.  Joints  were  measured 
in  a  road  outcrop  of  the  Eldon  formation  which  consists  of  grey,  medium  to  coarse 
crystalline  dolomites  showing  conspicuous  mottling  and  pseudo-algal  structures.  Beds 
are  from  1  foot  to  2  1/2  feet  thick.  The  area  is  located  on  the  east  limb  of  the  Bow 
River  anticline. 
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GEOMETRIC  ANALYSIS 


Introduction 

In  order  to  facilitate  the  classification  and  correlation  of  the  mesoscopic 
structures,  three  orthogonal  fabric  axes  a,  b,  and  c  were  defined  in  each  area. 
Rules  adopted  in  the  selection  of  the  fabric  axes  were  as  follows  (Turner  and  Weiss, 
1963,  p.  88): 

a)  In  folded  beds  the  fold  axis  =  b;  b  is  normal  to  the  plane  ac;  a  is  horizontal 

(Fig.  1). 

b)  In  uniformly  dipping  beds  the  plane  of  symmetry  =  ac;  the  normal  to  ac  lying  in 
the  bedding  plane  =  b;  a  is  horizontal. 

The  relationship  of  the  fabric  axes  to  uniformly  dipping  beds  is  taken  to  be  identical 
to  that  of  the  fabric  axes  to  the  limbs  of  folds  (Fig.  1).  The  choice  of  fabric  axes 
is  completely  free  of  any  kinematic  or  dynamic  implications. 


Bow  River  Area 

Folding.  Six  anticlines  in  the  competent  sandstone  units  of  the  Cardium  formation 
were  examined  in  detail  (Fig.  6).  The  folds  are  bilaterally  symmetric  across  the 
axial  plane.  Kananaskis  anticline  and  Horseshoe  anticline  approach  orthorhombic 
symmetry,  while  all  other  folds  display  a  definite  monoclinic  symmetry. 

Kananaskis  anticline  (Fig.  9),  Horseshoe  anticline  (Fig.  10),  and  Flattop 
anticline  (Fig.  13)  are  open  folds  with  axial  planes  that  are  steeply  inclined  to  the 
southwest.  Oldfort  anticline  (Fig.  11),  Cutoff  anticl ine  (Fig.  12),  and  Tight  anti¬ 
cline  (Fig.  14)  have  thrust-faulted  northeastern  limbs  which  display  steep  dips.  The 
east  limb  of  Tight  anticline  is  overturned  to  the  southwest.  Axial  planes  of  these 
folds  show  a  definite  southwesterly  dip. 

The  bedding  poles  of  each  fold  were  plotted  on  a  stereographic  projection 
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(Fig.  9-14).  For  each  fold  the  bedding  poles  tend  to  fall  along  one  great  circle  on 
the  stereonet.  Cylindrical  folding  requires  that  all  bedding  poles  lie  in  the  plane 
normal  to  the  fold  axis.  In  the  folds  analyzed,  the  poles  tend  to  scatter  about  the 
ideal  great  circle.  The  scatter  may  represent  the  deviation  of  the  folding  from  truly 
cylindrical  folding.  Large  scale  irregular  warping  of  the  bedding  surfaces,  as 
observed  at  several  localities,  would  tend  to  confirm  this  assumption.  On  the  other 
hand,  cross-stratification  in  the  sandstone  units  may  hve  influenced  bedding  plane 
determinations  sufficiently  to  account  for  the  scatter. 

All  the  folds  can  be  considered  to  be  essentially  planar  cylindrical  folds. 
There  appears  to  be  no  apparent  thickening  of  the  sandstone  units  in  the  axial  regions, 
or  thinning  in  the  limbs.  The  folds  are  of  the  concentric  type. 

-intersections  were  found  to  yield  spurious  results  when  used  for  the 
determination  of  the  fold  axis.  The  results  obtained  varied  widely  and  showed  a 
strong  dependence  on  the  positions  of  the  measurements  on  the  folded  surface,  and 
on  the  number  of  bedding  attitudes  measured  on  each  limb.  When  several  thousand 
intersections  were  used,  results  usually  corresponded  closely  with  value  for  the  {2  -axis 
obtained  from  point  diagrams. 

The  attitudes  of  the  fold  axes  (B)  were  determined  from  the  stereographic 
plots  of  the  bedding  planes  (Fig.  9-14).  The  trends  of  anticlinal  fold  axes  range 
from  N  26°  W  to  N  40°  W.  All  folds  plunge  to  the  northwest  at  angles  from  7  to  20 
degrees.  The  plunges  of  the  folds  increase  progressively  from  west  to  east. 

Bedding  -plane  si  ickensides_.  Slickensides  along  bedding  planes  throughout  the  area 
trend  from  N  42°  E  to  N  80°  E.  The  average  trend  of  the  slickensides  is  perpendicular 
to  the  mean  fold  axis  of  the  area.  At  some  localities  more  than  one  generation  of 
slickensides  can  be  differentiated  on  the  same  bedding  plane;  their  directions  usually 
do  not  differ  by  more  than  15  degrees. 
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Bedding  plane  si ickensides  represent  a  stress  release  along  bedding  planes 
during  folding.  The  presence  of  several  generations  of  movement  implies  discontinuous 
strain  during  folding. 

Joints  have  not  been  displaced  by  the  bedding-plane  slip. 

Joints.  The  maximum  number  of  joint  sets  is  observed  in  the  most  westerly  domain  of 
the  area  investigated.  Four  sets  of  joints  can  be  differentiated  in  the  west  limb  of 
Kananaskis  antic!  Ine  (Fig.  15).  Because  the  joints  in  this  domain  are  well  exposed 
and  typical  of  the  joints  found  throughout  the  area,  they  will  be  discussed  in  detail. 

a)  Two  conjugate  hkl  joint  sets  (N  32  E/81  SE  and  N  77  W/82  NE)  have  a  dihedral 
angle  of  72  degrees  and  an  acute  bisectrix  trending  N  68°  E.  The  N  32  E/81  SE 
set  shows  abundant  evidence  of  shear-type  movement  and  the  fractures  are  smooth 
and  planar.  The  N  77  W/82  NE  set  shows  only  rare  indications  of  shear-type 
movement  and  is  the  most  common. 

b)  An  hOI  joint  set  (N  30  W/82  NE)  consists  of  irregular  fractures  which  are 
nevertheless  large  and  prominent.  Evidence  of  normal-type  movement  is  abundant. 
In  rare  cases  fractures  belonging  to  this  set  show  plumose  markings  on  the  fracture 
surface. 

c)  A  Okl  joint  set  (N  66E/88  SE)  contains  some  planar  fractures  up  to  50  feet  long. 
Evidence  of  movement  along  the  joints  is  rare  and  can  be  both  of  the  shear-type 
and  the  normal -type.  Plumose  markings  on  the  joint  surface  can  be  observed.  The 
joint  set  parallels  the  bisectrix  of  the  conjugate  hkl  sets. 

Only  three  of  the  fracture  sets  observed  on  the  west  limb  of  the  anticline  are  present 
on  the  eastern  limb  (Fig.  16).  The  Okl  set  is  found  to  be  absent.  The  dihedral  angle 
between  the  conjugate  hkl  joint  sets  has  decreased  to  60  degrees. 

Two  domains  on  the  western  limb  of  Florseshoe  anticline  illustrate  the 
influence  of  lithology,  bedding  thickness,  and  outcrop  conditions  on  the  fracture 


70 


N 


1 80  ° 


FIG-  15 


N 


.  V  |  • 


71 


development.  The  domains  are  located  in  structurally  equivalent  positions.  One 
domain  includes  fractures  measured  in  the  upper  and  middle  sandstone  units  (Fig.  17) 
of  the  Cardium  formation,  while  the  other  represents  fractures  in  the  lower  sandstone 
unit  (Fig.  18).  The  mean  orientation  of  a  very  strong  maximum,  due  to  one  of  the 
hkl  sets  (N  83  W/87  NE),  is  essentially  identical  in  both  domains.  Marked  differences 
are  apparent  in  the  weaker  maxima.  The  complementary  hkl  set  is  nearly  absent  in 
the  lower  sandstone  unit.  The  strike  of  bedding  can  be  seen  to  change  by  90  degrees 
across  the  crest  of  Florseshoe  anticline  (Fig.  19).  The  orientations  of  the  joint  sets 
change  only  slightly  across  the  same  interval.  The  strike  of  the  bedding  has  no 
influence  on  the  orientation  of  the  joints  (Fig.  17-21). 

The  reproducibility  of  results  was  tested  on  the  eastern  limb  of  Horseshoe 
anticline.  Joints  were  measured  in  equivalent  structural  positions  and  stratigraphic 
units  on  both  banks  of  the  Bow  River  (Fig.  20  and  21).  Circles  of  confidence  about 
the  mean  orientations  of  the  two  domains  intersect  for  each  joint  set.  The  two 
orientation  diagrams  are  therefore  equivalent. 

Three  joint  sets  are  present  in  all  domains  on  Horseshoe  anticline.  The 
dihedral  angle  between  conjugate  hkl  sets  decreases  from  48  degrees  in  the  west  limb 
to  33  degrees  in  the  east  limb. 

In  the  west  limb  of  Oldfort  anticline  the  upper  and  middle  sandstone  units 
show  different  fracture  patterns.  Two  conjugate  hkl  joint  sets  in  the  upper  sandstone 
unit  (Fig.  23)  are  being  represented  by  a  single  Okl  joint  set  in  the  middle  sandstone 
unit  (Fig.  24).  The  conjugate  sets  intersect  with  a  dihedral  angle  of  33  degrees. 

Only  two  maxima  can  be  differentiated  in  the  western  limb  of  Cutoff  anit- 
cline  (Fig.  25  and  26)  and  the  east  and  west  limbs  of  Flattop  anticline  (Fig.  27  and  28). 
The  maximum,  occupying  the  position  of  the  two  conjugate  hkl  sets  found  west  of 
these  domains,  shows  an  unusually  large  dispersion  in  strike.  The  maximum  is  inter¬ 
preted  to  represent  two  merging  maxima.  Two  conjugate  hkl  sets  intersecting  at 
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low  dihedral  angles  are  thought  to  be  present  in  these  localities. 

Joints  on  the  west  limb  of  Tight  anticline  fall  into  an  hOI  and  a  Okl  set 
(Fig.  29  and  30). 

Joint  intersections.  Intersections  of  the  two  hkl  joint  sets,  and  intersections  of  the 
orthogonal  hOI  and  Okl  joint  sets  were  calculated  for  selected  domains.  The  inter¬ 
sections  were  found  to  form  maxima  about  the  normal  to  bedding.  Intersections  of 
conjugate  hkl  joint  sets  show  less  dispersion  than  intersections  of  the  Okl  and  hOI  sets. 

Variation  in  dihedral  angle.  A  plot  of  the  dihedral  angle  between  conjugate  hkl  sets 
versus  the  distance  perpendicular  to  the  reguional  strike  reveals  a  linear  relationship 
(Fig.  31).  For  dihedral  angles  smaller  than  35  degrees  the  maxima  of  the  two  hkl 
sets  merge,  and  the  determination  of  the  mean  fracture  planes  and  the  dihedral 
angle  between  them  becomes  impossible. 

Slickensides  on  joint  planes.  Nearly  15  per  cent  of  the  observed  joint  planes  show 
signs  of  movement.  All  conjugate  hkl  joints  and  Okl  joints  show  shear-type  slicken¬ 
sides.  These  slickensides  tend  to  parallel  the  intersection  between  the  joints  and  the 
bedding  plane.  No  conclusions  could  be  reached  about  the  directions  of  movement, 
because  the  sense  of  movement  deduced  from  the  roughness  of  the  slickensides  was  often 
contradictory  for  two  parallel,  adjacent  joints.  hOI  joints  invariably  show  normal- 
type  slickensides. 


Red  Deer  River  Area 

Two  orthogonal  joint  sets  predominate  in  the  Red  Deer  anticline.  The  mean 
fracture  planes  of  an  hOI  joint  set  strike  between  N  45°  W  and  N  64°  W.  Dips  of 
the  joints  are  highly  dependent  on  bedding  plane  dips  (Fig.  32-37).  Mean  fracture 
planes  of  a  Okl  set  strike  between  N  22°  E  and  N  36°  E.  A  third  joint  set,  which 
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is  only  weakly  developed,  can  be  seen  m  both  the  upper  and  lower  sandstone  units 
where  beds  dip  over  30  degrees  (Fig.  32  and  35).  The  mean  fracture  plane  of  this 
set  trends  N  50°  E. 

Okl  maxima  are  found  to  be  much  stronger  than  hOI  maxima  in  all  domains 
on  Red  Deer  anticline.  The  predominance  of  the  Okl  joints  can  be  at  least  partially 
attributed  to  outcrop  conditions.  Exposures  in  the  area  fall  along  a  line  roughly 
perpendicular  to  the  fold  axis,  and  Okl  joints  are  much  better  exposed  than  hOI  joints. 
Okl  joints  show  abundant  evidence  of  shear-type  movement.  The  sense  of  movement 
is  indeterminate,  hOI  joints  rarely  display  si ickensiding,  but  movement  along  them  is 
of  the  normal  type. 


Cripple  Creek  Area 

Joints  in  the  Cripple  Creek  area  are  found  to  be  irregular  and  frequently 
curved.  The  poles  of  the  joints  define  a  well  developed  girdle  with  three  weak 
maxima  (Fig.  38).  Two  conjugate  hkl  joint  sets  (N  85  E/80  NE  and  N  35  E/75  SW) 
exhibit  predominantly  shear-type  movement.  An  hOI  joint  set  (n  44  W/48  NE)  shows 
normal  movement  which  appears  to  have  taken  place  preferentially  in  a  plane  striking 
N  20°  W  (Fig.  39). 


South  Creek  Area 

Joints  were  measured  at  two  locations  in  equivalent  structural  positions 
on  the  ast  limb  of  a  large  syncline.  Three  sets  of  joints  are  present  in  the  western 
domain,  while  only  two  joint  sets  are  found  to  the  east.  A  Okl  joint  set  (N  40° 
E/88  SE)  and  an  hOI  joint  set  (N  33°  W/72  NE)  are  common  to  both  localities 
(Fig.  40).  Shear-type  slickensides  on  the  Okl  joints  indicates  sinistral  relative 
movement  and  normal -type  movement  has  occurred  on  the  hOI  joints  (Fig.  41). 

An  hkl  joint  set  (N  77°  W/72  NE)  observed  in  the  western  domain  appears  to  bear 
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no  relationship  to  the  regional  structures. 

Exshaw  Area 

Th  ree  very  conspicuous  joint  sets  occur  in  the  massive  beds  of  limestones. 
Two  conjugate  hkl  joint  sets  (N  37°  E/70  SE  and  N  73°  E/81  NW)  form  very  sharply 
defined  maxima  and  intersect  at  a  dihedral  angle  of  46  degrees.  An  hOI  joint  set 
(N  43°  W/51  NE)  with  a  larger  dispersion  than  the  hkl  sets  is  also  present  (Fig.  42). 
Notable  at  this  locality  Is  the  complete  absence  of  si  Ickensiding  and  the  presence  of 
zones,  1-2  feet  thick,  in  which  the  fracture  frequency  is  considerably  higher  than  in 
the  surrounding  rock  body. 


The  Gap  Area 

The  joints  In  the  limestones  of  the  Banff  formation  show  large  dispersions 
in  strike,  but  two  major  joint  sets  can  be  differentiated  (Fig.  43).  A  Ok  I  joint  set 
(N  48  E/81  SE)  with  shear-type  si  Ickensiding  makes  an  angle  of  60  degrees  with  a 
set  of  joints  (N  13  W/  62  NE)  which  show  no  indications  of  movement.  Curved 
tension  gashes  filled  with  crystalline  calcite  are  found  to  trend  parallel  to  bedding 
and  indicate  that  shearing  has  taken  place  along  bedding.  Joints  cut  these 
structures  without  being  displaced  by  the  bedding  slip, 

Roche  MSette  Area 

Although  fractures  In  the  Pal  User  limestones  at  this  locality  are  abundant 
and  regular,  they  are  found  to  have  only  poor  preferred  orientations.  A  Ok  I  joint 
set  (N  44  E/  78  NE)  is  definitely  present,  but  no  other  joint  sets  can  be  distinguished. 
The  normal  to  bedding  forms  the  axis  to  a  wide  zone  Into  which  the  remainder  of  the 
joints  fa  1 1 . 
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Morro  Peak  Area 

The  degree  of  preferred  orientation  of  joints  both  in  the  Pal  User  and  the 
Alexo  formations  is  low  in  this  area.  A  Ok  I  joint  set  can  be  defined,  but  the 
validity  of  other  concentrations  is  dubious  (Fig.  44  &  45).  Fractures  in  the 
Palliser  formation  fall  into  a  well  defined  girdle. 

Mount  Murchison  Area 

Joints  in  the  Eldon  dolomites  on  the  east  limb  of  the  Bow  River  anti¬ 
cline  fall  into  two  distinct  maxima  (Fig.  46).  The  mean  fracture  planes  are 
orientated  N  43°  W/  80°SW  and  N  60°  E/  89°  SE.  Shear-type  slickensides 
are  present  on  both  joint  sets.  (Fig.  47), 

Summary 

The  geometry  of  the  fractures  in  the  Foothills  of  the  Canadian  Rocky 
Mountains  bears  a  definite  relationship  to  the  geometry  of  the  folds  and  faults. 
The  observed  joints  usually  belong  to  an  hOI  set,  okl  set,  or  pair  of  hkl  sets  of 
fracture  planes.  The  Okl  and  hkl  fractures  are  found  to  be  generally  smooth, 
planar,  and  extensive.  Slickensiding  on  their  surfaces  indicates  that  shear- 
type  movement  has  taken  place  along  them.  An  hOI  joint  set  is  found  to  be 
present  in  all  domains  investigated  in  the  Foothills.  The  hOI  fractures 
usually  are  rough  and  discontinuous  in  appearance  and  show  large  deviations 
in  their  strike  directions.  Slickensiding  along  the  hOI  surfaces  invariably 
indicates  normal-type  movement. 

Investigations  in  the  Eastern  and  Main  Ranges  show  that  the 
relationship  between  the  macroscopic  structure  and  the  joints  is  essentially 
the  same  as  in  the  Foothills.  An  exception  are  joints  in  the  Eastern  Ranges 
east  of  Jasper  which  have  unusually  large  dispersions. 
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Joints  formed  at  a  post-folding  stage  in  the  tectonic  history  of  the 
Rocky  Mountains.  None  of  the  joints  are  found  to  be  displaced  by  bedding 
plane  slip  or  have  been  deformed  by  folding.  Attempts  to  establish  the 
relative  ages  of  different  joint  sets  by  their  cutting  relationships  in  the  out¬ 
crop  were  unsuccessful.  Joints  of  every  set  are  found  to  cut  joints  of  every 
other  set.  Indirectly  the  age  relationships  can  be  inferred  from  the  dispersions 
shown  by  the  different  joint  sets.  First-formed  joints  should  exhibit  lower 
dispersions  than  those  formed  after  the  contiguity  of  the  rock  body  was  lost 
due  to  fracturing,  hkl  and  Okl  joint  sets  usually  show  less  dispersion  than 
the  hOI  joint  sets.  The  dispersion  of  lines  of  intersection  of  conjugate  hkl 
sets  is  also  markedly  lower  than  for  intersections  of  Okl  and  hOI  sets. 

The  results  of  a  study  of  the  dispersion  shown  by  joint  sets  were 
negative.  Dispersion  is  not  the  direct  result  of  lithology,  bedding  thickness, 
structural  position,  outcrop  conditions,  or  intensity  of  deformation. 

Fracture  frequencies  are  dependent  on  the  bedding  thickness  and 
joints  show  a  closer  spacing  in  thin  beds.  No  differences  in  the  fracture 
frequencies  of  the  limbs  and  crests  of  folds  could  be  detected.  Fracture 
frequency  does  not  appear  to  be  a  function  of  the  magnitude  of  the  tectonic 
forces. 
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KINEMATIC  AND  DYNAMIC  INTERPRETATIONS 

Introduction 

The  geometric  analysis  of  the  joints  revealed  the  existence  of  a  very 
consistent  fracture  pattern  in  the  Foothills  of  the  Canadian  Rocky  Mountains. 
Investigations  in  the  Eastern  and  Main  Ranges  indicate  that  a  similar  pattern  may 
also  exist  in  those  regions.  The  mechanisms  of  formation  of  joints  in  all  the  areas 
investigated,  as  well  as  the  stress  field,  must  therefore  have  been  similar. 

Any  theory  explaining  the  formation  of  the  joints  has  to  account  for  the 
following  observations: 

a)  the  geometry  of  the  joint  sets; 

b)  the  movements  on  the  joint  surfaces; 

c)  the  changes  in  dihedral  angles  between  conjugate  joint  sets; 

d)  the  variability  in  the  number  of  joint  sets  present; 

e)  the  development  of  the  joints  at  a  late  stage  in  the  tectonic  history. 

Principal  Stress  Trajectories 

The  most  persistent  geometric  property  of  the  joints  investigated  in  each 
domain  is  their  tendency  to  fall  along  a  girdle,  which  can  be  approximated  by  a 
great  circle  on  the  stereographic  net.  The  presence  of  a  girdle  indicates  that  the 
joints  are  cozonal  about  an  axis,  which  is  given  by  the  girdle  axis.  In  cases  where 
the  girdle  axis  coincides  with  the  normal  to  the  bedding,  the  joints  are  perpendicular 
to  the  bedding.  In  the  majority  of  areas  investigated  the  girdle  axes  are  found  to 
plunge  approximately  2-10  degrees  steeper  than  the  normal  to  the  bedding.  Two 
areas  in  the  Eastern  Ranges,  along  the  Bow  River,  are  found  to  be  exceptional  in 
having  girdle  axes  inclined  at  a  lower  angle  than  the  normal  to  bedding  (Fig.  42 
&  Fig.  43). 
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Photoelastic  experiments  recently  carried  out  by  Bell  and  Currie  (1964, 
pp.  33-51)  showed  that  the  maximum  principal  stress  trajectories  are  subparallel  to 
the  upper  and  lower  boundaries  of  competent  units.  In  one  of  their  experiments  a 
single  competent  plastic  member  was  surrounded  by  incompetent  gelatine  and  folded 
concentrically  under  uniaxial  compression.  The  maximum  principal  stress  trajectories 
were  found  to  run  within  the  competent  member.  In  the  limbs  of  the  fold  the 
maximum  principal  stress  trajectories  were  found  to  be  inclined  at  a  slightly  lower 
angle  than  the  bedding  (Bell  and  Currie,  1964,  p.  43,  Fig.  5b). 

The  experimental  model  described  above  is  almost  completely  analogous 
to  the  situation  observed  in  the  Cardium  formation  of  the  Foothills.  Single  competent 
sandstone  members  enclosed  in  incompetent  shales  have  been  folded  concentrically. 
Principal  stress  trajectories  subparallel  to  the  boundaries  of  the  competent  sandstone 
members  explain  the  formation  of  fractures  virtually  perpendicular  to  the  bedding. 
Furthermore,  principal  stress  trajectories  plunging  at  slightly  lower  angles  than  the 
bedding  dip  would  result  in  joints  cozonal  about  an  axis  plunging  at  a  higher  angle 
than  the  normal  to  bedding. 

Similar  conclusions  have  been  reached  by  other  workers.  Charlesworth  and 
Evans  (1962,  pp.  358-359)  studied  the  deviation  of  slaty  cleavage  from  the  mean 
plane  and  demonstrated  that  the  maximum  principal  stress  axis  acted  parallel  to  the 
competent  members  in  the  stratigraphic  succession.  Hancock  (1964,  p.  184)  investi¬ 
gated  the  late-formed  joints  in  the  folds  of  South  Pembrokeshire  and  concluded  that 

the  trace  of  the  <5~  axis  must  have  paralleled  bedding  during  the  formation  of  the 

X 

joints. 

The  occurrence  of  girdle  axes  inclined  at  a  lower  angle  than  the  normal 
to  bedding  is  confined  to  two  localities  in  the  Eastern  Ranges  (Fig.  42  &  Fig.  43). 
Both  stations  are  located  in  thick  competent  units  that  have  been  faulted  along 
steeply  dipping  thrust  planes.  The  deflection  of  the  maximum  principal  stress 
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trajectories  toward  the  thrust  plane,  rather  than  parallel  to  the  boundaries  of  the 
competent  unit,  would  account  for  the  geometry  of  the  joints. 

Significance  of  Variations  in  Dihedral  Angle 

Experimental  studies  have  shown  that  the  dihedral  angle  between  conjugate 
fracture  sets  will  increase  with  increasing  confining  pressure  (Karman,  1911). 
Muehlberger  (1961,  pp.  211-219)  suggested  a  modification  to  the  Mohr  envelope 
of  fracture  to  account  for  this  variability  of  dihedral  angle.  Furthermore,  he 
theorized  that  over  a  critical  stress  interval  a  gradational  change  from  a  single 
extension  fracture  to  two  conjugate  shear  fractures  should  occur.  Duschatko  (1953) 
has  been  the  only  author  to  encounter  such  a  relationship  in  the  field. 

In  the  present  study  the  exact  situation  postulated  by  Muehlberger  (ibid.) 
is  found  to  exist,  and  the  change  of  the  dihedral  angle  is  found  to  vary  directly  with 
distance  across  the  structural  trend  (Fig.  31).  Two  conjugate  joint  sets  at  high 
dihedral  angle  in  the  western  portion  of  the  Bow  River  area  are  seen  to  have  decreasing 
dihedral  angles  to  the  east.  Eventually  the  two  sets  are  found  to  be  replaced  by  a 
single  joint  set. 

Using  the  Mohr  envelope  and  the  modifications  suggested  by  Muehlberger 
(ibid.),  the  observations  made  in  the  Bow  River  area  can  be  explained.  At  the  time 
of  joint  formation  the  confining  pressures  were  highest  in  the  western  portions  of  the 
area  and  large  differential  stresses  existed.  Conjugate  shear  joint  sets  with  high 
dihedral  angles  formed  (Fig.  48,  point  A).  The  confining  pressure,  and  as  a  result 
also  the  differential  stresses,  decreased  to  the  east.  The  radius  of  curvature  of  the 
circle  began  to  approach  that  of  the  parabolic  vertex  of  the  envelope.  Although 
stress  circles  were  still  tangent  to  the  envelope  at  two  points,  the  dihedral  angles 
progressively  decreased  with  decreasing  differential  stresses  (Fig.  48,  point  B). 

Even  farther  east,  the  stress  circle  eventually  had  the  same  radius  of  curvature  as 
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the  parabola  forming  the  vertex  of  the  envelope.  Only  one  set  fractures  formed 
at  that  point  (Fig.  48,  point  C). 

The  confining  pressure  is  a  function  of  the  weight  of  overburden.  The 
western  portion  of  the  Bow  River  area  must  at  one  time  have  been  under  a  higher 
overburden  than  areas  to  the  east.  Overburden  in  the  Red  Deer  River  and  South 
Creek  areas  must  have  been  equivalent  to,  or  less  than,  that  in  the  eastern  portions 
of  the  Bow  River  area.  At  present,  all  these  locations  are  at  approximately  the  same 
topographic  elevations. 

The  trace  of  the  McConnell  thrust  fault  is  found  to  be  3  miles  west  of  the 
western  border  of  the  Bow  River  area.  The  attitude  of  the  thrust  plane  can  be  seen 
to  change  from  steeply  dipping  in  the  west  to  nearly  flat  in  the  vicinity  of  Mount 
Yamnuska,  located  a  few  miles  north  of  Kananaskis  Falls.  At  one  time  the  hanging 
wall  of  the  McConnell  thrust  must  have  extended  much  farther  east  than  it  does  today. 
The  rocks  in  the  Bow  River  area  formed  part  of  the  footwall  of  the  McConnell  thrust. 
Somewhere  in  the  region  east  of  its  present  trace,  the  thrust  must  have  cut  up  in  the 
stratigraphic  section  in  order  to  die  out  eventually.  In  that  case,  the  thrust  sheet 
must  have  thinned  to  the  east.  Rocks  of  the  same  structural  horizon  in  the  footwall 
would,  therefore,  have  been  overlain  by  a  thcker  thrust  sheet  than  those  in  the  east. 
This  model  would  account  for  the  apparent  differences  in  confining  pressure.  If  the 
thrust  plane  had  cut  up  in  the  stratigraphic  section  at  a  constant  angle,  then  the 
linear  relationship  observed  between  the  dihedral  angles  and  the  distance  across  the 
regional  strike  would  be  explained. 

This  hypothesis  can  be  tested  by  applying  the  results  derived  in  the  Bow  River 
area  to  some  of  the  other  areas  investigated.  The  Red  Deer  River  and  South  Creek 
areas  are  located  12  and  8  miles  east  of  the  McConnell  thrust  respectively.  At  those 
distances  the  angles  between  conjugate  shear  sets  would  be  either  so  low  that  the  sets 
could  not  be  differentiated,  or  only  one  extension  set  would  form.  In  both  areas  that 
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is  the  observed  situation.  The  Cripple  Creek  area  lies  approximately  5  miles  east 
of  the  McConnell  thrust,  and  two  conjugate  shear  sets  with  dihedral  angle  of  47 
degrees  should  form  in  this  area.  Two  weak  shear  joint  sets  are  observed  at  dihedral 
angle  of  49  degrees  .  An  approximate  correlation  between  the  confining  pressure 
under  which  the  joints  formed  and  the  distance  from  the  present  trace  of  the  McConnell 
thrust  seems  therefore  possible.  Only  further  observations  could  substantiate  this 
finding. 


Origin  of  Joints 

During  the  main  tectonic  phase  in  the  Rocky  Mountains  all  principal 
stresses  were  compressional,  and  the  formation  of  extensional  joints  is  inconceivable. 
The  development  of  shear  joints  is  dependent  on  two  conditions: 

a)  the  intermediate  stress  axis  has  to  be  vertical; 

b)  the  differential  stresses  have  to  be  large  enough  to  cause  shearing  in  the  rock. 

Price  (1959,  pp.  153-155)  demonstrated  that  both  of  these  conditions  cannot  be 
satisfied  at  the  same  time  during  the  main  phase  of  compression.  Geological  evidence 
yields  further  proof  for  the  late  formation  of  the  joints. 

At  the  end  of  the  period  of  thrusting  and  folding  the  following  stresses  acted 

in  the  areas  investigated: 

acting  from  northeast  to  southwest  = 

I  x 

<^2  acting  from  southeast  to  northwest  =  C 

^2  acting  vertical  to  bedding  =  c5^ 

where  ,  <5*  and  6  0  are  the  maximum,  intermediate,  and  minimum  principal 
stresses  respectively.  C  ^  acts  in  the  ac  fabric  plane,  C5  ^  in  the  be  fabric  plane,  and 
®  is  perpendicular  to  the  ab  fabric  plane 

Not  all  the  stresses  had  been  dissipated  during  the  thrusting  and  folding 
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because  the  rocks  have  a  fundamental  strength.  The  residual  stresses  remained 
stored  in  the  rocks  and  acted  parallel  to  the  boundaries  of  the  competent  units 
(Fig.  49,  point  A). 

The  two  lateral  stresses  ^  and  Cf~  consisted  of  two  components: 

x  y 

a)  stresses  due  to  gravitational  loading  equal  to 


1 


m  -  1 


b)  tectonic  stresses  which  are  C  and 

a 


* 

C~  z 


1  Cm  the  x  and  y  directions 
— -a  ' 

m 


respectively. 

m  is  equal  to  Poisson's  number. 

As  the  area  underwent  erosion,  the  amount  of  gravitational  loading  £T ^ 
decreased  and  resulted  in  a  proportional  decrease  in  the  gravitational  components  of 
the  lateral  stresses.  In  addition,  the  rock  body  underwent  lateral  strain  in  response 
to  uplift.  Lateral  tensile  stresses  equal  to  ahlf  the  change  in  load  developed  (Price, 
1959,  p.  158).  As  a  result,  the  lateral  stresses  decreased  at  a  more  rapid  rate  than 
during  unloading.  At  some  point  (Fig.  49,  point  B)  the  vertical  stress  changed 
to  the  intermediate  principal  stress  and  Cf^  became  the  minimum  principal  stress. 
One  of  the  conditions  necessary  for  the  formation  of  shear  joints  was  now  satisfied. 
With  further  uplift,  the  ratio  of  3  increased  rapidly.  When  this  ratio  reached 

a  critical  value,  two  conjugate  hkl  shear  joint  sets  developed  (Fig.  49,  point  C). 

In  areas  of  low  confining  pressures,  conjugate  sets  of  low  dihedral  angle  or  only  a 
Okl  extension  set  formed  (Fig.  49,  point  D). 

With  the  formation  of  the  shear  sets,  or  the  single  set  of  extension  joints, 
a  large  fraction  of  the  residual  stresses  was  dissipated.  The  maximum  principal  stress 
was  now<5~z,  the  intermediate  principal  stress  <5^,  and  the  minimum  principal 
stress  CS*  .  Continued  erosion  and  uplift  resulted  in  tensile  stresses  that  eventually 
decreased  Cf  to  the  point  where  it  became  tensile.  Shortly  thereafter,  a  set  of 
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hO!  extension  fractures  formed  perpendicular  to  the  minimum  principal  stress  axis. 

The  formation  of  the  hOI  joints  immediately  released  the  tensile  stresses;  <5^  became 
the  minimum  principal  stress. 

Further  uplift  resulted  in  tensile  stresses  in  the  a-direction.  In  areas  where 
a  Ok  I  fracture  set  or  two  conjugate  joint  sets  at  small  dihedral  angle  had  already 
formed,  these  stresses  resulted  in  the  opening  up  of  the  previously  formed  joints. 

In  regions  of  large  dihedral  angles,  a  Okl  fracture  set  would  form  in  the  absence  of  a 
suitable  plane  along  which  the  tensile  stresses  could  dissipate.  This  is  observed  to 
be  the  case  in  the  west  limb  of  Kananaskis  anticline  where  two  conjugate  hkl  joint 
sets,  a  hOI  joint  set,  as  well  as  a  Okl  joint  set  are  present. 

The  movements  observed  on  the  joint  surfaces  supports  this  hypothesis. 

Shear  movements  are  observed  on  all  the  Okl  joints  and  conjugate  sets  of  hkl  joints, 
illustrating  their  close  genetic  relationship.  Normal-type  movement  prevades  on 
the  hOI  joint  sets  and  agrees  with  an  extensional  origin. 
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Direction  cosines  of  counting  locations  used  on  Program  913107-001  . 
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APPENDIX  11 


Domains  of  the  Bow  River  (Bl-20),  Red  Deer  River  (RD1-7),  Cripple 
Creek  (Cl),  and  South  Creek  (SI -2)  areas.  The  structural  position  of  each 
domain  is  described  where  possible  with  reference  to  folds,  and  the  strati¬ 
graphic  position  in  terms  of  the  lower  (L),  middle  (M)  and  upper  (U)  sandstone 
units  (ss)  of  the  Cardium  Formation.  The  number  of  measurements  (N),  mean 
strike  (S),  mean  dip  (D),  radius  of  the  95  percent  confidence  circle  about 
the  mean  (CR),  and  precision  parameter  (K)  for  each  joint  set  in  each  domain 
are  given;  for  bedding  only  the  mean  strike,  mean  dip  and  number  of  measure¬ 
ments  are  shown . 
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Domain 

Location 

Structure 

N 

S 

D 

C.R. 

K 

B1 

West  limb 

bedding 

29 

175 

14  SW 

Kananaski's  Anticline 

hoi 

79 

150 

82  NE 

2.7 

36 

okl 

68 

66 

88  SE 

2.0 

71 

L  M  U  ss 

hkl 

66 

32 

81  SE 

1.9 

89 

hkl 

80 

104 

82  NE 

4.4 

14 

B2 

East  limb 

bedding 

13 

107 

11  NE 

Kananaski's  Anticline 

hoi 

52 

153 

87  SW 

3.9 

26 

hkl 

100 

41 

87  SE 

1.8 

61 

L  M  ss 

hkl 

85 

101 

86  SW 

1.4 

122 

B3 

West  limb 

bedding 

24 

4 

16  NW 

Horseshoe  Anticline 

hoi 

87 

160 

79  NE 

3.0 

27 

hkl 

98 

50 

83  SE 

2.9 

26 

L  M  U  ss 

hkl 

54 

97 

87  NE 

3.9 

16 

B4 

West  limb 

bedding 

43 

7 

17  NW 

Horseshoe  Anticline 

hoi 

92 

174 

79  NE 

3.6 

17 

hkl 

30 

35 

75  SE 

7.9 

12 

L  M  ss 

hkl 

126 

100 

85  NE 

1.7 

67 

B5 

Crest 

bedding 

45 

51 

11  NW 

Horseshoe  Anticline 

hoi 

55 

165 

89  NE 

2.8 

47 

hkl 

75 

52 

80  SE 

3.0 

31 

L  ss 

hkl 

101 

91 

86  SW 

1.7 

66 

B6 

East  limb 

bedding 

37 

96 

15  NE 

Horseshoe  Anticline 

hoi 

103 

163 

85  SW 

2.1 

46 

hkl 

70 

63 

80  SE 

2.3 

54 

L  M  ss 

hkl 

65 

99 

78  SW 

1.6 

122 

B7 

East  limb 

bedding 

18 

no 

14  NE 

Horseshoe  Anticline 

hoi 

93 

162 

86  SW 

2.6 

33 

hkl 

39 

51 

82  SE 

3.1 

14 

L  M  U  ss 

hkl 

113 

88 

81  SE 

1.8 

53 

B8 

East  limb 

bedding 

25 

95 

15  NE 

Horseshoe  Anticline 

hoi 

80 

165 

87  SW 

2.3 

48 

hkl 

64 

56 

84  SE 

2.4 

57 

L  M  U  ss 

hkl 

92 

92 

79  SW 

1.7 

73 

B9 

West  1  imb 

bedding 

20 

36 

12  NW 

Oldfort  Anticline 

hoi 

102 

163 

86  NE 

2.4 

35 

hkl 

82 

59 

83  SE 

2.2 

52 

U  ss 

hkl 

45 

91 

85  SW 

2.5 

74 

BIO 

West  1  imb 

bedding 

40 

166 

38  SW 

Oldfort  Anticline 

hoi 

104 

153 

60  NE 

3.3 

18 

M  ss 

okl 

122 

62 

82  SE 

2.4 

18 

116 


Bll 

West  limb 

bedding 

22 

33 

14  NW 

Cutoff  Anticl ine 

hoi 

120 

155 

83  NE 

2.0 

44 

M  ss 

ok! 

113 

67 

83  SE 

2.5 

28 

B12 

West  limb 

bedding 

11 

41 

13  NW 

Cutoff  Anticline 

hoi 

68 

152 

88  NE 

2.2 

61 

M  ss 

okl 

117 

68 

84  SE 

1.4 

92 

B13 

West  limb 

bedding 

20 

32 

14  NW 

Cutoff  Anticl  ine 

hoi 

58 

147 

81  NE 

4.  1 

22 

hkl 

46 

20 

77  SE 

5.8 

14 

L  ss 

hkl 

140 

83 

87  SE 

2.3 

27 

B14 

West  limb 

bedding 

11 

173 

22  SW 

Flattop  Anticline 

hoi 

90 

160 

70  NE 

2.3 

42 

U  ss 

okl 

127 

61 

80  SE 

2.7 

22 

B15 

West  limb 

bedding 

14 

8 

20  NW 

Flattop  Anticline 

hoi 

88 

169 

76  NE 

3.0 

26 

L  M  ss 

okl 

133 

68 

83  SE 

3.0 

18 

B16 

Crest 

bedding 

10 

40 

10  NW 

Flattop  Anticline 

hoi 

71 

153 

86  NE 

2.5 

45 

L  ss 

okl 

168 

60 

81  SE 

2.9 

16 

B17 

East  limb 

bedding 

8 

71 

19  NW 

Flattop  Anticline 

hoi 

102 

144 

88  SW 

2.6 

30 

L  ss 

okl 

131 

55 

80  SE 

3.3 

15 

B18 

East  1  imb 

bedding 

7 

76 

16  NW 

Flattop  Anticline 

hoi 

122 

150 

89  NE 

2.7 

24 

L  ss 

okl 

120 

57 

77  SE 

3.2 

17 

B19 

West  1  imb 

bedding 

10 

32 

22  NW 

Tight  Anticl  ine 

hoi 

74 

153 

80  NE 

2.0 

72 

L  ss 

okl 

92 

57 

78  SE 

1.5 

99 

B20 

West  1  imb 

bedding 

9 

163 

48  SW 

Tight  Anticl  ine 

hoi 

75 

154 

59  NE 

5.7 

9 

L  ss 

okl 

111 

52 

80  SE 

3.3 

17 

RD1 

Crest 

bedding 

19 

50 

6  NW 

hoi 

70 

125 

87  SW 

2.8 

38 

L  ss 

okl 

137 

38 

85  SE 

2.3 

29 

RD2 

Crest 

bedding 

24 

34 

7  NW 

hoi 

116 

128 

89  SW 

2.0 

44 

U  ss 

okl 

171 

35 

85  SE 

1.4 

59 

RD3 

East  limb 

bedding 

39 

118 

14  NE 

hoi 

115 

124 

73  SW 

3.8 

13 

L  ss 

okl 

155 

22 

89  SE 

1.5 

54 

/ 

117 

RD4 

East  limb 

bedding 

hoi 

L  ss 

okl 

RD5 

East  limb 

beddinq 

hoi 

U  ss 

okl 

RD6 

East  limb 

bedding 

hoi 

L  ss 

hkl 

hkl 

RD7 

East  limb 

bedding 

hoi 

U  ss 

okl 

Cl 

L  U  ss 

bedding 

hoi 

okl 

SI 

L  ss 

bedding 

hoi 

hkl 

hkl 

47 

143 

34  NE 

119 

131 

56  SW 

3.2 

17 

151 

26 

74  NW 

2.5 

30 

57 

138 

29  NE 

102 

127 

63  SW 

2.0 

48 

176 

30 

82  NW 

1.6 

44 

15 

149 

41  NE 

2.3 

89 

136 

55  SW 

3.7 

17 

139 

19 

63  NW 

2.5 

35 

43 

52 

79  NW 

10.1 

5 

32 

151 

38  NE 

2.1 

70 

117 

63  SW 

4.5 

15 

154 

20 

61  NW 

2.3 

24 

12 

147 

38  SW 

103 

136 

48  NE 

3.1 

21 

148 

59 

88  NW 

3.7 

11 

13 

134 

16  SW 

117 

147 

72  NE 

1.3 

101 

127 

40 

88  SE 

1.8 

48 

33 

103 

72  NE 

2.8 

81 

